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0> ■ Nucleosynthesis in ONeMg novae has been investigated with the wide ranges of three 

^ . parameters, i.e., the white dwarf mass, the envelope mass at ignition, and the initial com- 

position. A quasi-analytic one-zone approach is used with an up-to-date nuclear reaction 
O ' network. The nucleosynthesis results show correlation with the peak temperatures or 

^ • the cooling timescales during outbursts. Among the combinations of white dwarf and 

ci ■ envelope masses which give the same peak temperature, the explosion is more violent 

for a lower white dwarf mass owing to its smaller gravitational potential. Comparison of 
the nucleosynthesis results with observations implies that at least two-third of the white 
dwarf masses for the observed ONeMg novae are ~ I.IMq, which are significantly lower 
than estimated by previous hydrodynamic studies but consistent with the observations of 
V1974 Cyg. Moreover, the envelope masses derived from the comparison are ^ 10~^Mq, 
which are in good agreement with the ejecta masses estimated from observations but 
significantly higher than in previous hydrodynamic studies. With such a low mass white 
dwarf and a high mass envelope, the nova can produce interesting amounts of 7-ray 
emitters ''Be, ^^Na, and ^^Al. We suggest that V1974 Cyg has produced ^^Na as high as 
the upper limit derived from the COMPTEL survey. In addition, a non-negligible part 
of the Galactic ^^Al may originate from ONeMg novae, if not the major contributors. 
Both the future INTEGRAL survey for these 7-ray emitters and abundance estimates 
derived from ultraviolet, optical, and near infrared spectroscopies will impose a severe 
constraint on the current nova models. 

Subject headings: novae, cataclysmic variables — nuclear reactions, nucleosynthesis, 
abundances — white dwarfs 
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1. INTRODUCTION 

A classical nova has been thought to be a ther- 
monuclear runaway of hydrogen-rich gas accumu 
l ated onto a w hite dwarf in a close binary system 
Gehrz et 



1982 



al. 



1998 and references 
thc pin). Recent observations show that about 30% 
of well-studied events are classified as oxygen-neon- 



it is still significantly lower than observational esti- 
mates. Third, there has been no consensus on the 
mixing mechanism between the white dwarf matter 
and the accreted gas, though a few hypotheses such as 
diffusion, shear mixing, and convective overshooting 
have been proposed (Prialnik & Kovetz 1984; Kuttei 



fc Sparks 1987; [ben, Fujimoto, fc MacDonald 1991 



magnesium (ONcMg) novae. Observationally, ONcMg' 



Glasner, Livne, & Truran 1997; Kcrcek, Hillebrandt 



novae are characterized by strong line emissions in 
neon and other intermediate-mass elements like mag- 
nesium, aluminum, silicon, and sulfur in their ejected 



shetTs~( Livio & I'ruran lyy4 ). I he presence of these 



elements implies that the accumulated gases must 
have been substantially enriched through the dredge- 
up [from the ONeMg cores. 



ONcMg novae have been suggested to be a promis- 
ing production site of 7-ray emitters ^Be, ^^Na, and 



26 Al (Btarrfield, Truran, & Sparks 1978; 


Neiss & Tru- 


ran 199C; Nofar, Shaviv, & Starrfield 1991 


Starrfield 


3tal. 1993; Coc et al. 1995; Politano et al. 1995; Her- 


nanz et al. 199fi|; Wanajo et al. 1997a, b; 


Jose, Her- 


tianz, & Coc 1997; Jose & Hernanz 1998; 


Starrfield et 


il. 19981). However, the following three uncertainties 



confront us when studying nucleosynthesis in ONeMg 
novae. First, the mass of the ONeMg white dwarf 
is not constrained from theoretical models any more 



than ~ 1.1 - 1.4MtT,, which results from the 8 - lOM, 



stellar evolution models (Nomoto 1984, 1987; [ben 
'k Tutukov 1985 ). On the other hand, only a few 
observational estimates of the white dwarf masses 



sported ( 


Paresce et al. 1995 




Krautter 


Retter, Leibowitz, & Ofek 1997 


). Sec- 



ond, there is a serious disagreement on the accreted 
masses onto white dwarfs between observational esti- 
mates and current theories. The ONeMg white dwarfs 
in previous hydrodynamic studies accumulate a few 
10~^Mq of the envelope masses at ignition (Politano 



3t al. 1995; 3tarrficld et al. 1998; Jose & Her 



aanz 1998). On the other hand, the estimated ejecta 



masses of QU Vul, V838 Her, and V1974 Cyg are 



~ io-^-io-^Mq ( 


Taylor et al. 1987; 


Greenhouse et 


il. 1988; Baizar et al. 1992; Woodward et al. 1992|; 


Pavelin et al. 1993 


; Shore et al. 1993 


; Baizar et al. 


1996 


; Vanlandingham et al. 1996|; Woodward et al. 


1997 


), which are 10—100 times larger than theoretical 



estimates. Starrfield et al. (1998) have shown that the 
envelope mass increases with decreasing mass accre- 



tion rate and wh i te dwarf luminosity (see a lso Prialnik 
^ Kovetz 1995[ [Kovetz fc Prialnik 1997| ). However, 



& Truran 1998a, b). Furthermore, the metallicity es- 
timates for the observed ejecta of ONeMg novae show 
a wide spread between 0.09 and 0.86 in ma s s fraction 
(iLivio fc Truran 1994| ; [Pohtano et al. 1995| ; Starrfield 



et al. 1 



The initial composition of an envelope 
may significantly affect the nucleosynthesis result as 



well as the energetics of the outb urst ( [Kovetz fc Pri- 
alnik 1997|; [Jose fc Hernanz 1998|). 



The purpose of this study is to examine nucle- 
osynthesis in ONeMg novae with the wide ranges of 
three parameters: the white dwarf mass, the enve- 
lope mass, and the mixing ratio of the core-surface 
matter into the envelope. In § H, we describe our 
quasi-analytic nova models and an updated nuclear 
reaction network. We then, in § compare the nu- 
cleosynthesis results for one sequence with a previ- 
ous hydrodynamic calculation. In § ^, we constrain 
the ranges of white dwarf and envelope masses, com- 
paring the nucleosynthesis results with observational 
abundance estimates; where the effect of changing the 
initial composition is considered. Finally, the 7-ray 
line emissions from ^Be, ^^Na, and ^^Al are discussed 
in § 

2. METHOD OF CALCULATION 
2.1. Nova Model 

Our nova models are based on the quasi-analytic 
approach for the hydrogen shell fiash on a white dwarf 



(Sugimoto fc Fujimoto 1978; Fujimoto 1982a, b). The 



temperature and density structures of an envelope are 
obtained analytically for a given set of a white dwarf 
mass (Mwd) and an envelope mass (Mgnv), on the 
assumption that the spherical envelope expands in 
hydrostatic equilibrium. We have constructed mod- 
els for 49 sets of Mwd (1.05 - 1.35Mq) and Menv 
(10~6 — 1O~^M0). The former corresponds to the 
masses of ONeMg cores which results from 8 — IOMq 
stellar evolu tions (Nomoto 1984, 1987; |lben fc Tu" 



tukov 1985), and the latter covers those both from 



theories 10 ^ 




- lO-^M©; [Truran et al. 1977 




Poli- 


tano et al. 1995 




Starrfield et al. 1998|; Jose fc Her- 
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nanz 1998) and from observations 10 '^Mq; Tay- 


lor 
al. 


et al. 1987; Greenhouse et al. 198 


3; ^aizar et 


1992|; [Woodward et al. 1992; Pavehn et al. 1993|; 


shore et al. 1993 


; Saizar et al. 1996; 


Vanlanding- 


hail et al. 1996]; 


Woodward et al. 1997 


). The dots 



cal calculations are performed, while squares, trian- 
gles, and stars are taken from hydrodynamic studies 
by Politano et al. (1995, hereafter PSTWS95), Star- 
rfield et al. (1998, hereafter STWS98), and Jose & 
Hernanz (1998, hereafter JH98). The solid lines show 
the mass accretion rates onto the white dwarfs re- 
quired for each set of (Mwd, Afcnv), calculated by 
Fujimoto (1982b). These are in reasonable agree- 
ment with those by PSTWS95, STWS98, and JH98 
(■~ 10~^° — IO^^'Mq yr~^), but somewhat overesti- 
mated since the luminosities of white dwarfs are ne- 
glected and the radii are assumed to be of Chan- 
drasekhar (see Figure ||) in Fujimoto (1982b). Note 
that no outburst is achievable by an accreting white 
dwarf below the dashed line due to the high accretion 
rate. It is obvious that a rather low accretion rate (or 
a low luminosity of the white dwarf) is required to 
obtain a massive envelope such as ~ 10^^ — 10~^Mq 
as expected by observations. 

The quasi-analytic nova model has been elaborated 
by Sugimoto & Fujimoto (1978) and Fujimoto (1982a, 
b). Let us discuss the model in some detail, since it 
can characterize the nova burst very well. The pres- 
sure and the density at the base of the envelope are 
expressed in terms of Mwd and Mem '■ 

nv 



where i?wD is the radius of the white dwarf, is a 
homologous invariant defined by 

d\nP _ GMp 
d In r rP 

Hereafter the subscript 'b' denotes a quantity at the 
base of the envelope. The flatness parameter / in the 
equations (|l|) and (^ decreases monotonically as the 
shell flash proceeds: 



3-N 



where 



iV+ 1 
V 



(0 < a; < 1) . 



(3) 



(4) 



The value of /b denotes the degree of the 'flatness' of 
the envelope. For /b ~ 1 {xh ~ 0), the envelope is 
thin and strongly degenerated, and thus is flat. On 
the other hand, for fh^O (a;b ^ 1), the envelope is 
thick and non degenerate, and thus is spherical. The 
polytropic index TV in the equations (^ and (^) is 
defined by 

N _ dlnp 
N+1 " dlnP ' 
and (p, q) is the incomplete beta function defined 

by 



Bx (p, q) 



tP-l (1 _ dt (0 < X < 1) . 



N is assumed to be adiabatic and constant through- 
out the envelope, but vary with time. The effect of the 
spatial variation in N is quite small for a typical con- 
vective envelope ( Fujimoto 1982a ). The value of N 
is approximately 1.5 at the beginning of a shell flash, 
and approaches 3 at the end due to the increasing 
radiation pressure. 

The shefi flash starts with /b 1 (a^b ^ 0). The 
envelope is then heated up by nuclear burning to a 
thermal runaway, and cools down when /b decreases 
to ^ (a;b 1). The equations (|l|) and (||) are valid 
if 

Ur 

(5) 



1 



•C 1 



is satisfled, where 



U = 



din A// 47rr^p 



dlnr M 
is another homologous invariant. This condition is 
violated only near the last phase of the shell flash 
(/b ~ 0). At this phase, major nuclear reactions are 
frozen out except for the pp-chain, the CNO cycle, 
and /3"'"-decay. Thus, our nucleosynthesis results may 
not be signiflcantly affected. 

Figure]^ illustrates contours for Phi jh and ph/Vbfh 
in the M-wD-Menv space. These are the the proper 
quantities for each set of (Mwd, Menv)- The stronger 
dependence of the former on Mwd is due to the higher 
power of i?wD as seen in the equation (|^) and (||). 
The temperature at the base of the envelope Tb can 
be calculated by solving the equation of state with the 
use of equations (|l|) and (||). The spatial variations 
of the pressure, the density, and the temperature are 
given when the condition (||) is satisfied, by 

P (x) = Pb ■ ■ ' 



Xh 



1 - Xb 



3 



p{x) 
T{x) 



= Ph\ — 



= Th — 



1 — X 

1 - 

{N+1)V 



-N 



1 



l-Xh 



-{N+1)V 



where V = d\nT/d\nP is assumed to be adiabatic 
and constant throughout the envelope, but vary with 
time (on the deviation from constant V, see Fujimoto 
1982a). The value of x decreases monotonically with 
increasing radius, approaching zero at the surface of 
the envelope. The surface radius R is given when the 
condition (h) is satisfied, by 



R = 



1 - Xb' 



(6) 



Now we know the envelope structure completely. 

The progress of a shell flash is derived by energy 
conservation, 

dt (T) ' 

where En is the nuclear energy generation rate per 
unit mass, s is the specific entropy which is spatially 
constant in the convective envelope, and (T) is the 
mass averaged temperature over the envelope. The 
energy inflow from the white dwarf and loss from 
the photosphere are neglected, being much smaller 
than the nuclear energy during the explosive hydro- 
gen burning. The time variation of Xh is then calcu- 
lated from the equations (|l|), (||), and (0) with the 
use of the equation of state. The expansion velocity 
of the envelope Uoxp is derived from the equation (^) 
as 

R dxh 
'^^^ 1 — x\, dt 

Each calculation is started with the initial tempera- 
ture Tb = 5 X lO'' K, and ceased when the nuclear lu- 
minosity decreases to the Eddington luminosity where 
no further heavy elements are synthesized. 

The MwD~^WD relation is derived for an isother- 
mal core (2 x 10^ K) consists of oxygen, neon (= 
5:3), and partially degenerate electron gase s in- 



clu ding the effect of th e Coulomb interaction ( [chi 



u fc Kitamura 1994 ), as shown in Figure ||. The 



solid line denotes our results and the triangles are 
taken from PSTWS95 and STWS98. Our results 
are between those of carbon and magnesium white 
dwarfs by Hamada & Salpeter (1961), and somewhat 
smaller than by PSTWS95 and STWS98. A varia- 
tion of i?wD significantly influences the density due 



to pb c>c i?wD~'^ as seen in the equation (||), much 
more than the temperature (oc i?wD~^)- Note that 
the ONe white dwarf is unable to increase its mass 
beyond 1.38M0 because the electron capture on ^"Ne 
and ^^Mg triggers the collapse (denoted by a dot on 
the solid line; Nomoto 1984, 1987). 

2.2. Nuclear reaction network and initial com- 
position 

The nuclear reaction network used in this work 
contains 87 stable and proton-rich isotopes from hy- 
drogen to calcium (Table |l|), including all relevant 
nuclear reactions and weak interactions. The reac- 
tion *B(p, 7)^C, which can be a sink for the ^Be pro- 
duction (Boffin et al. 1993), is also included. The 
ground and isomeric states of ^^Al take longer than 
the mean lifetime of the i somer (::± 9.2 s) to be equi- 



librated for ^ 4 X 10^ K ([Ward fc Fowler 1980| ). The 
peak temperatures in the models responsible for the 
observed ONeMg novae may be less than 4 x 10^ K 
as will be discussed in § ||. Thus, the two states are 
separated as different isotopes. The nuclear reaction 
rates are taken from Thielemann et al. (1995). They 
are based on the rates by Caughlam & F owler (1988) 



those calculated by a statistical model ( Fruran et al 



1987), and the latest experimental data (Van Wormei 
et al. 1994, etc.). We also include new reaction rates 
by Herndl et al. (1995) and Ihadis et al. (1996). 



The rate ^^Si(p, 7)^^? (Herndl et al. 1995) may have 



a special importance, being 10"^ — 10^ times larger 
than the previous one in the typical nova tempera- 
ture range. The rates ^^Mg(p, 7)^^A1 and ^^Al(p, 
7)^^Si ( Iliadis et al. 1996| ) may be also of importance 



for °A1 production, though the latter involves a large 
uncertainty. In our computations, all nuclear reaction 
rates are mass-averaged over the envelope except for 
/J^-decay which does not depend on density and tem- 
perature. 

The initial composition of an envelope is assumed 
to be a mixture of the solar composition gas and the 
dredged-up matter from the surface of the ONeMg 
white dwarf. The solar abun dances are adopted from 
Anders & Grevesse (1989), and the abundances of 



the ONeMg core m atter from plashimoto, Iwamoto, 
fc Nomoto 
(Table g). 



(1993)1 for the 1.35 Mq ONeMg core 
As can be seen in Table § O : Ne : 



Mg « 10 : 6 : 1, which is in good agreement with 
those in Nomoto and Hashimoto (1988) for Mwd = 
1.26, 1.36M0 and Ritossa, Garcia, & Iben (1996) for 
MwD = 1.2M0. This implies that the composition 



4 



of an ONeMg core does not significantly depend on 
its mass. The mass fraction of the dredge-up matter 



from the ONeMg core in the envelope Xwd, which 



and density profiles taken from hydrodynamic studies 
( Starrfield, Truran, fc Sparks 197^ ; ^tarrfield, Sparks " 



is the third parameter in this study, is of importance 
on the nucleosynthesis results as will be discussed in 
However, 



abundance estimates in the obser- 
vations of nova ejecta involve large uncertainties as 
pointed out by Livio & Truran (1994). The estimated 
metallicities of the six observed ONeMg nova ejecta 
range widely (see Table |^) and, unfortunately, differ- 
ent authors have pr ovided different values even for 
the identical events (IWiUiams et al. 1985|; |Snijders| 



Starrficld 1994; Austin et al. 199e; 3aizar et al 



et il. 1987 ;_ Saizar et al. 1992|; Andrea, Drechsel 



19Si 6|; [Vanlandingham et al. 1996| [Vanlandingham 
Staprfield, fc Shore 1997). In addition, no consensus 



has been achieved in theoretical modeling how and 
when the core- matter mixes into the envelope ( Pri 



aln k & Kovctz 1984; Iben, Fujimoto, & MacDon 
aid 1991|; [Kuttcr fc Sparks 1987|; [Glasner, Livne, fc 
Tr 4ran 1997| ; [Kercck, Hillebrandt, fc Truran 19984 
b). Thus, we examine all the combinations of (Mwd, 
Monv) for XwD = 0.1 (case A), 0.4 (case B), and 0.8 



(case C), which cover observational uncertainties in 
abundance determinations. I'hc initial compositions 
for each case are given in Table 2.2. 



Comparison with nucleosynthesis by a hy- 
drodynamic model 

Up to now, a number of works on nucleosynthesis 



in ONeMg novae have been performed ( 


Hillebrandt 


fc Thielemann 1982; 


Weiss fc Truran 1990 


; Nofar, 


Bhaviv, fc Starrficld 1991 


, and references therein). 



Their nova models were, however, based on one-zone 
envelopes, using the spatially constant temperature 



Table 1: Nuclear Reaction Network Employed 



Element 


4 ■ 


A 

^max 


Element 


A ■ 

^inin 


A 

^inax 


H 


1 


2 


Na 


20 


23 


He 


3 


4 


Mg . . . . 


21 


26 


Li 


7 


7 


Al 


22 


27 


Be 


7 


7 


Si 


24 


30 


B 


8 


11 


P 


27 


31 


C 


9 


13 


S 


28 


34 


N 


13 


15 


CI 


31 


37 





14 


18 


Ar 


32 


38 


F 


17 


19 


K 


35 


39 


Ne 


18 


22 


Ca 


36 


40 



Shaviv 198^). Coc et al. (1995) have studied ^^Na 



and ^Al production in ONeMg novae with another 
semi-analytic method (MacDonald 1983). Their nova 
model and ours give similar envelope structures in 
temperature and density. However, our model in- 
cludes the effect of the partially degenerate and rela- 
tivistic electron gas, while Coc et al. (1995) treated 
electrons as the ideal gas. The electron degeneracy 
can not be neglected in the early phase of outbursts. 
Hernanz et al. (1996) and Jose, Hcrnanz, fc Coc 
(1997) have also examined nucleosynthesis in novae 
with the use of a hydrodynamic method. However, 
they focused on ^Li or ^^Al production, and gave only 
a few synthesized isotopes in their papers. 

Hence, we compare our model with sequence 6 in 
STWS98 to see the differences of nucleosynthesis be- 
tween the quasi-analytic and hydrodynamic methods. 
The nova model in STWS98 was identical to that of 
PSTWS95, except that the former included the up- 
dated nuclear reacti on rates ( Van Wormer et al. 1994 



Herndl et al. 1995 ) and OPAL opacity tables ( Igle- 
sias fc Rogers 19931 ). In addition, STWS98 employed 
a lower white dwarf luminosity and a lower mass ac- 
cretion rate to obtain a more massive ignition enve- 
lope. Furthermore, an important change was that 
STWS98 used a longer mixing length of (2 — 3) x pres- 
sure scale height. We do not compare our results with 
JH98 which has studied nucleosynthesis in ONeMg 
(and CO) novae using a hydrodynamic code, since 
the white dwarf radii are not presented. Their re- 
sults showed, however, similar trends to PSTWS95 
and STWS98. We use the same initial composition, 
the nuclear reaction rates, Mwd {— 1.25Mq), M^nv 
(= 4.5 X lO-^M©), and i?wD as STWS98 for com- 
parison. Note that the nucleosynthesis results in this 
work are obtained for the whole envelope, while those 



Table 2: Abundances of the ONeMg Core at the Sur- 
face 



Nucleus 


Mass Fraction 


Nucleus 


Mass Fraction 




3.95E-02 


24Mg 


4.20E-02 


160 


5.42E-01 


25Mg 


6.29E-03 


20Ne 


3.31E-01 


26Mg 


4.57E-03 


21 Ne 


2.87E-03 


27A1 


1.25E-02 


22Ne 


1.34E-Q3 


28Si 


2.46E-03 


23Na 


1.65E-02 
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Table 3: Initial Compositions of the Envelope by Mass 



X^^u 0.1 0.4 0.8 



p 


6, 


,36E- 


-01 


4 


.24E- 


-01 


1, 


.41E- 


■01 


D 


4, 


.33E- 


-05 


2 


.88E- 


-05 


9, 


,62E- 


-06 


3He 


2, 


.64E- 


■05 


1 


.76E- 


■05 


5, 


.87E- 


■06 


^Hc 


2, 


,48E- 


-01 


1 


.65E- 


-01 


5, 


.51E- 


-02 


^Li 


8, 


.43E- 


-09 


5, 


.62E- 


-09 


1, 


.87E- 


-09 




4, 


.26E- 


■09 


2 


.84E- 


■09 


9, 


.46E- 


■10 


12c 


6, 


,69E- 


-03 


1 


.76E- 


-02 


3, 


.22E- 


■02 


13c 


3, 


,29E- 


-05 


2 


.19E- 


-05 


7, 


,31E- 


-06 




9, 


,96E- 


■04 


6 


.64E- 


■04 


2. 


.21E- 


■04 


15N 


3, 


,93E- 


-06 


2 


.62E- 


-06 


8, 


.74E- 


■07 


160 


6. 


.28E- 


-02 


2 


.22E- 


-01 


4, 


.35E- 


-01 


1^0 


3, 


,50E- 


■06 


2 


.34E- 


■06 


7, 


.79E- 


■07 


180 


1, 


.95E- 


-05 


1 


.30E- 


-05 


4, 


.34E- 


■06 


19p 


3, 


.65E- 


■07 


2, 


.43E- 


-07 


8, 


,11E- 


-08 


20Ne 


3. 


.45E- 


■02 


1 


.33E- 


■01 


2, 


.65E- 


■01 


2iNe 


2, 


,90E- 


-04 


1 


.15E- 


-03 


2, 


.29E- 


-03 


22Ne 


2, 


,51E- 


■04 


6 


.15E- 


-04 


1, 


,10E- 


-03 


23Na 


1, 


.68E- 


■03 


6 


.60E- 


■03 


1, 


.32E- 


■02 


24Mg 


4, 


.66E- 


-03 


1 


.71E- 


-02 


3 


.37E- 


■02 


25Mg 


6, 


.89E- 


-04 


2 


.55E- 


-03 


5, 


,04E- 


-03 


26Mg 


5, 


,27E- 


■04 


1 


.88E- 


■03 


3 


.67E- 


■03 


27 Al 


1, 


.31E- 


-03 


5 


.05E- 


-03 


1, 


.OOE- 


■02 


28Si 


8, 


.34E- 


-04 


1 


.38E- 


-03 


2, 


,10E- 


-03 


29Si 


3, 


.09E- 


■05 


2 


.06E- 


■05 


6, 


.86E- 


■06 


30Si 


2, 


.12E- 


-05 


1 


.41E- 


-05 


4, 


.71E- 


■06 


31p 


7, 


.35E- 


-06 


4, 


.90E- 


-06 


1, 


.63E- 


-06 


32s 


3, 


.57E- 


■04 


2 


.38E- 


■04 


7, 


.93E- 


■05 




2, 


,90E- 


-06 


1 


.94E- 


-06 


6 


.45E- 


-07 


34g 


1, 


.68E- 


-05 


1 


.12E- 


-05 


3, 


,74E- 


-06 


35C1 


2, 


,28E- 


■06 


1 


.52E- 


■06 


5. 


.07E- 


■07 


37C1 


7, 


.70E- 


-07 


5 


.13E- 


-07 


1, 


.71E- 


-07 


36 Ar 


6, 


.98E- 


-05 


4, 


.65E- 


-05 


1, 


.55E- 


-05 


38 Ar 


1, 


,39E- 


■05 


9 


.24E- 


■06 


3 


.08E- 


■06 


39k 


3, 


,13E- 


-06 


2 


.08E- 


-06 


6 


.95E- 


-07 


40Ca 


5. 


.40E- 


-05 


3 


.60E- 


-05 


1, 


,20E- 


-05 
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in STWS98 for the ejected matter. Thus, STWS98 
may strongly reflect the composition of the outer re- 
gion. Figure ^ shows the ratios of isotopes (dots) 
and elements (triangles) between ours and STWS98 
(sequence 6). Our calculation obtains a higher peak 
temperature (~ 3.17 x 10^ K) than that of STWS98 
(~ 3.00 X 10^ K), since the latter model ignited hydro- 
gen one zone above the base so that the envelope is 
effectively thinner (see STWS98). The prominent un- 
derproduction of several isotopes hke ^^N, ^^O, ^^Ne, 
22Na (and perhaps ^^Na not shown in STWS98; see 
PSTWS95 for instance), ^'^Mg, and ^^Mg is due to 
our assumption of a fully convective one-zone enve- 
lope. Since these isotopes are rather fragile against 
the (p, 7) or (p, a) reactions, they decrease signifi- 
cantly even at the late phase of the outburst. In con- 
trast, these isotopes were able to survive in STWS98, 
escaping from the hotter convective region into the 
cooler radiative region at the late phase. Especially, 
-'^^N is extremely fragile against the (p, a) reaction, 
thus being underproduced by more than 5 orders of 
magnitude in this work. As a result, nitrogen (mostly 
^^N in this work) is also underproduced compared to 
STWS98 in which ^^N is dominant. On the other 
hand, carbon (^^C and ^'^C) are significantly overpro- 
duced, transferred from ^^N. We should be careful on 
these differences in comparing the nucleosynthesis re- 
sults with observations. However, both results are in 
excellent agreement for other isotopes, and especially 
for elements (except for carbon and nitrogen) which 
are more important for comparison with observations. 

4. Nucleosynthesis in ONeMg novae 
4.1. Nuclear flovirs in the N—Z plane 

In this section, we present some important aspects 
of nucleosynthesis in ONeMg novae, referring to the 
results of several (Mwd, -^^cnv) models. Figure ^ 
shows the final abundances and the net nuclear flows 
in the N-Z plane. The size of a circle denotes the 
mole fraction of the isotope defined by 1^ = Xi/Ai 
in the logarithmic scale. The initial composition is 
shown by dotted circles. The net nuclear fiow of a 
reaction from the i-th to j-th isotope, defined as 



Y^{i^3)-Yj{3^i) 



dt . 



metallicity of the ejecta estimated from observations 
(see 'Z' in Table |). 

Figure ^ shows the peak temperature at the base 
Tpcak, the cooling timescale r defined as the dura- 
tion from the peak to its half in temperature, the 
peak nuclear energy generation rate per unit mass 
£pcak, and the ejection velocity Wej in the Af-wD~A?onv 
space. Here, Wcj is defined as the expansion velocity 
Wexp when it equals the escape velocity Wcsc (for the 
models denoted by circles). For the models denoted 
by crosses in which Uexp is below Vcac throughout the 
calculations, Uej is replaced with Uoxp at the maxi- 
mum. As seen in Figure |^, r has a weaker depen- 
dence on Mwd than Tpeak, while the trend of £pcak 
is similar to Tpeak- As a result, among the models 
of the same peak temperature, the explosion is more 
violent for the smaller Mwd due to its smaller gravi- 
tational potential. This is also seen in the panel of Woj , 
which shows the similar trend to r in the 7\fwD~-Menv 
space. In order to obtain the fast ejection velocities 
such as > 1000 km s~^ as derived by recent observa- 



is denoted by the length of an arrow in the same scale. 
The mixing ratio ATwd is assumed to be 0.4 (case B) 
throughout this section, which is close to the average 



tions ( |Gehrz ct al. 1998 and references therein), the 
cooling timescale must be <^ 1000 s where the /?"'"- 
decay of I'^O (r ~ 102 s) and (r ~ 176 s) plays 
an important role. 

4.1-1. Low temperature sequences 

For the model (Mwd/M©, Mc„v/Mq) = (1.10, 
lO^'''^), the initially present ^''Mg is entirely trans- 
ferred to silicon, even though Tpeak is as low as 

2 X 10* K (Figure |). In contrast, the initial ^'^Ne 
remains mostly unburnt, though minor nuclear flows 
appear through the Ne-Na cycle. A part of the ini- 
tial is converted to ^^O, ^^C, i^C, and ^^N. The 
HCNO cycle is active near the peak in temperature, 
turning to the CNO cycle as the temperature de- 
creases. Thus, almost all ^^N is eventually converted 
to ^''N, i^c, and "C. Note that, for the models with 
Tpeak 2 X 10* K, Wexp is too Small to overcome Wesc 
as seen in Figure ^. 

4-1.2. Moderate temperature sequences 

The nucleosynthesis results for (Mwd/M©, Menv 
/Mq) = (1.15, 10-4-°) and (1.35, lO'^-^) (hereafter 
N1540B and N3555B, respectively) differ significantly, 
regardless of their mostly same Tpeak (— 2.9 x 10* K) 
as seen in Figure |. This can be explained as fol- 
lows. Figure ^ shows the time variations of Tb and 
e for each model. The cooling timescale for N1540B 



7 



(r ~ 190 s) is more than one order shorter than for 
N3555B (r ~ 2400 s). This is a consequence of the 
weaker gravitational potential for N1540B owing to its 
smaller Mwd and thus its larger i?wD (Figure ||) . In 
addition, the nuclear energy generation rate remains 
as high as ^ 10^** erg g~^ s~^ even after the envelope 
expands and the temperature decreases to ^ 10^ K, 
owing to the /J^-decay of ^''O, ^^O, and other unsta- 
ble nuclei. As a result, the expansion of the envelope 
is accelerated and then the temperature drops fairly 
quickly, even when its structure is returning to the 
static configuration. In contrast, for N3555B, almost 
all the short-lived f3~^ -unstable nuclei have decayed at 
the late phase. Hence, the temperature drops slowly 
with the decreasing nuclear energy generation rate. 
The patterns of the temperature decreases are, there- 
fore, not similar between these models. The criti- 
cal cooling timescale between the slow (N3555B) and 
fast (N1540B) expansion is t 1000 s. The cooling 
timescale for N1540B is comparable to the /?"'"-decay 
lifetime of ^^O (= 176 s). As a result, ^^N survives the 
following (p, a) reactions and significantly enhances. 
For similar reasons, "'^^O, ^^Mg, and ^^Al are promi- 
nent in N1540B, while they are absent in N3555B. 
Note that the somewhat higher epoak in N1540B is 
due to the higher density at the base (Figure ||). 

It is noteworthy that the net nuclear flows of 
^'*Mg(p, 7)^^A1 have overcome the initial abundance 
of 24Mg for both N1540B and N3555B (Figure |), ow- 
ing to substantial nuclear flux from the Ne-Na region. 
It implies that the initial amount of ^''Mg does not sig- 
nificantly affect the production of isotopes A> 24: for 
the models Tpoak ^ 3 x 10® K. Note that N1540B also 
obtains the significantly higher ejection velocity (~ 
2100 km s-i) than N3555B (~ 1200 km s^^). As seen 
in Figure |6[ for all the models with Tpoak ^ 3 x 10® K, 
VcKp exceeds Vcsc and obtains v^j > 1000 km s~^, 
which is in good agreement with recent observations 
of ONeMg novae. 

4-1.3. High temperature sequences 

For the models {Mwu/Mq, M^^v/Mq) = (1.20, 
10-* ") and (1.20, lO-^-^) (hereafter N2040B and 
N2035B, respectively), substantial nuclear fluxes ap- 
pear in the Ne-Na region because of their high Tpoak 
(~ 3.3 X 10® K and 4.2 x 10® K, respectively) as seen 
in Figure In addition, various nuclear paths open 
in the Mg-S region. The abundance of ^^Al is highly 
enhanced in N2040B due to the substantial nuclear 
flux from the Nc-Na region via ^^Na(p, 7)^''Mg. On 



the other hand, ^^Al is less abundant in N2035B be- 
cause of its higher peak temperature. Instead, ^®0, 
^^Na, and ^■^Na are highly enhanced in N2035B, since 
T (~ 9.5 s) is comparable to the /3'*'-decay lifetimes of 
i®Ne (2.4 s), 22Mg (5.6 s), and ^^Mg (16 s). For the 
extremely high temperature (Tpoak — 7.3 x 10® K) 
model (Mwd/Mq, M^^^/Mq) = (1.30, lO-^"), al- 
most all the initial ^"Ne is burned out, and the nuclear 
flow extends to calcium by the rp-process (Figure ^ . 
The leakage from the CNO cycle via the a-capture of 
^^O and ^^O appears, though its contribution to the 
heavy element production is negligible. 

4.2. Element and isotope production 

In this section, we discuss the global trends of ele- 
ment production and isotope ratios in the MwD-Afonv 
space, referring to the abundances of ONeMg nova 
ejecta estimated from recent observations. Table ^ 
shows the abun dances for the recent six ONeMg no- 



3el, & Starrfield 1994; 


Vanlandingham, Starrfield, & 


Shore 1997 


), V1370Aql(3nijdcrsctal. 1987; 


Andrea, 


Drcchscl, & Starrfield 1994), QU Vul (Saizar et al. 


1992; 


Andrea, Drechsel, & Starrfield 1994 


), V351 Pup 



(Saizar et al. 1996), V838 Her (Vanlandingham, Star- 



rfield, fc Shore 1997D , and V1974 Cyg ( [Austin ct al 
1996| ). Note that the abundances of the elements not 



presented in the above references are assumed to be 
zero, thus involving errors by a few percent. The av- 
erage metallicity for these ONeMg novae is ~ 0.43 by 
mass. The mixing ratio Xwd is, therefore, assumed 
to be 0.4 (case B) throughout this section. However, 
VI 370 Aql and V838 Her show significantly different 
metallicities from case B. The dependence on the ini- 
tial composition is discussed in § |4.3| . 

When temperature is higher than ^ 2 x 10® K, 
proton captures are fast enough to compete with the 
/3+-decay of various unstable isotopes. As a result, 
the nucleosynthesis results are significantly deviated 
from those in steady nuclear flows like the CNO and 
Ne-Na cycles. Figures |-0show the final abundances 
and isotope ratios by mass in the MwD--^env space. 
The abundances are shaded from white (0.1) to black 
(10"'°^) in the logarithmic scale (except for beryllium 
and boron) . In the rest of this paper, all abundances 
are given in mass fraction. As described below, we 
find that there exist two types of elements, namely, 
those correlated to Tpeak (e.g., oxygen, neon, and sul- 
fur) and to r (e.g., carbon, sodium, and magnesium). 
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Table 4: Observed ONeMg Nova Abundances 





H 


He 


C 


N 


O 


Nc 


Mg 


Al 


Si 


S 


Z 


V693 CrA^ 


2.8E-01 


3.2E-01 


5.1E-03 


8.4E-02 


1.2E-01 


1.7E-01 


7.6E-03 


3.4E-03 


2.6E-03 




4.0E-01 


V693 CrA^' 


1.6E-01 


1.8E-01 


7.9E-03 


1.4E-01 


2.1E-01 


2.7E-01 


1.8E-02 




6.9E-03 




6.5E-01 


V693 CrA'' 


3.9E-01 


2.0E-01 


4.3E-03 


8.0E-02 


7.5E-02 


2.3E-01 


2.9E-03 


1.9E-03 


8.7E-03 




4.1E-01 


V1370 Aql* 


4.9E-02 


8.8E-02 


3.5E-02 


1.4E-01 


5.1E-02 


5.2E-01 


6.8E-03 




1.8E-03 


l.OE-01 


8.6E-01 


V1370 Aql^ 


4.5E-02 


l.OE-01 


5.0E-02 


1.9E-01 


3.7E-02 


5.6E-01 


7.9E-03 




4.6E-03 




8.5E-01 


QU Vul= 


3.0E-01 


6.0E-01 


l.OE-03 


2.1E-02 


1.6E-02 


2.3E-02 


1.7E-03 




4.0E-02 




l.OE-01 


QU Vul^ 


3.3E-01 


2.7E-01 


9.6E-03 


7.4E-02 


1.8E-01 


8.7E-02 


3.7E-03 


9.9E-03 


3.2E-02 


1.2E-02 


4.0E-01 


V351 Pup'' 


3.8E-01 


2.4E-01 


5.9E-03 


7.4E-02 


1.9E-01 


l.lE-01 




4.3E-03 


1.9E-03 




3.8E-01 


V838 Hcr=' 


6.0E-01 


3.1E-01 


1.2E-02 


1.4E-02 


2.5E-03 


5.8E-02 








2.8E-03 


9.0E-02 


V1974 Cyg^ 


1.8E-01 


3.1E-01 


5.4E-02 


7.7E-02 


2.7E-01 


l.lE-01 










5.1E-01 



References: ^Williams et al. 1985, ^Andrea, Drcchscl, &: Starrficld 1994, ^Vanlandingham, Starrficld, & Shore 1997, ^Snijdor ct al. 1987, 
^Saizar ct al. 1992, "Saizar ct al. 1996, ''Austin ot al. 1969 



4-2.1. Beryllium and Boron 

As seen in Figure ||, the abundance of "^Be (in mass 
fraction) reaches ~ lO'^ for Tp^ak ~ 2.5 - 4 x 10^ K 
(Figure by the a-capture of the initially present 
^He. For the same Tpcak, the lower Af-yvD models pro- 
duce more ^Be than higher ones. This is due to the 
higher densities for the formers as seen in Figure ^. 
When density is less than ~ 10^ g cm^^ at temper- 
ature ~2 — 4xlO*K, the proto n capture of ^Be i s 
suppressed by its inverse reaction ( Bofhn ct al. 1993 ) . 
For Tpeak ^ 4 X 10^ K, ''Be decreases by its a-capture. 
As a result, the abundance of ^^B reaches ^ 10~^. For 
Tpcak ^ 6 X 10^ K, the abundance of ^^B decreases 
owing to the reaction ^^C{a, p). 

4.2.2. Carbon and nitrogen 

In the steady flow of the CNO cycle (<:, 2 x 10^ K), 



the l most abundant isotope is N and the isotope ra- 



tios arc determined by the nuclear reaction rates as 



12C/13C 
14N/15N 



A[i3c(p,7)] /A p2c(p,7)] ~2-4 
A[i5N(p,a)] /A ["N(p,7)] 
5000 - 50000. 



When temperature exceeds ~ 2 x 10* K, the CNO 
cycle is replaced with the HCNO cycle via ^^N(p, 
7)^^0(/3+i')^^N. The abundance patterns of the car- 
bon and nitrogen (Figure mainly depend on r (Fig- 
ure ^) as follows: (1) For r ^ 1000 s, the carbon 
and nitrogen isotopes show the typical feature of the 
steady CNO cycle, i.e., C/N < 1, ^^C/^^C ~ 3, and 
14^/15^ ^ 30000. (2) For t - 1000 s, however, these 
isotope ratios approach ^ 1, due to the /3+-decay life- 
times of i^N (~ 862 s) and ^^O (~ 176 s) compa- 
rable to the cooling timescale. The thermonuclear 
runaway ceases before most ^■^N (and some ^^O) de- 



cays, and thus the ratio C/N also reaches ^1. (3) For 
T <^ 1000 s, the thermonuclear runaway ceases dur- 
ing the active HCNO cycle where ^"^O and ^^O are 
abundant, resulting in C/N < 1. The ratio ^^C/^^C 
is unchanged (~ 3), while ^^N/-'^^N is significantly re- 
duced to ~ 0.1. 

The abundance of nitrogen is ~ 0.1 in the whole 
area of the MwD-M^nv space, regardless of the ra- 
tio ^^N/^^N ranging over 5 orders of magnitude. 
In contrast, the abundance of carbon ranges widely 
(~ 0.001 — 0.1) reaching its maximum at r ~ 1000 s, 
while the ratio ^^C/^^C is not significantly changed 
in the M^^-Mcnv space. The above results explain 
the abundance feature of the recent ONeMg novae 
(Table ^) , in which the abundance of carbon spreads 
widely 0.001 — 0.01) while that of nitrogen is ~ 0.1. 
Note that the abu ndance of nitrog en for QU Vul 
( Baizar et al. 1992 ) and V838 Her ( Vanlandingham 
et al. 1996 ) is as low as ~ 0.02, owing to the signifi- 
cantly lower metallicities (~ 0.1). For V838 Her and 
V1974 Cyg, the ratio of C/N is ~ 1 which is obtained 
by the models with r ^ 1000 s. 

It should be noted that our models may signifi- 
cantly underproduce ^^N, that causes the too large 
ratio C/N as discussed in § ^. This may be, how- 
ever, only the case in the models r ^ 1000 s. For 
the models r <^ 1000 s, the abundance of ^^N is not 
significantly reduced as described above, and thus the 
results may not be changed substantially. 

4.2.3. Oxygen and fluorine 

The abundance of oxygen is mainly correlated to 
Tpcak but is also dependent on r (Figure |o|), owing 
to the presence of three isotopes. The ratio ^^O/^^O 
has a clear correlation to Tpeak- It reaches the min- 
imum (~ 0.3) at Tpcak ~ 3 X 10* K, and is nearly 
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constant 3 for Tpoak <; 2 x 10^ K and - 10 
for Tpcak ^ 4 X 10* K), due to the different nu- 
clear reaction cycles (Figure ||). In contrast, the ratio 
16q^18q si^Q-^^s a, clear correlation with the cooling 
timescale (Figure |To| ), being significantly small for 
T ^ 100 s. As a result, the abundance of oxygen 
reaches ~ 0.03 - 0.1 for (1) Tpoak < 3 x 10* K {^^O 
and ^"^O are abundant) or for (2) r ^ 100 s (^®0 is 
abundant). Note that oxygen is always abundant in 
the models A/wd I.ISM©, where one of these con- 
ditions is satisfied. 

Fluorine (^^F) is not significantly enhanced in the 
all models (Figure The reason is that the reac- 
tion ^*F(p, 7)^^Ne, which is followed by the /J^-decay 
to ^^F, is much slower than ^^F{p, q;)^^0. The abun- 
dance of ^^F is ^ 10~'^ at most for t 10 s, which is 
comparable to the /J^-decay lifetime of ^^Ne (~ 25 s). 

The oxygen-rich ONeMg novae (~ 0.1 — 0.3 by 
mass) V693 CrA, QU Vul, V351 Pup, and V1974 Cyg 
(Table |) can be explained by the following models: 
(1) MwD < 1.15M0, (2) Tpeak < 2 X 10* K, or (3) 



r <j 10 s. On Llic oLhcr hand, VS3S Her is fairly 
geii poor (~ 6.6 X iu 
a rather massive model (Mwd 



IS iairly uxy- 
which could be explained by 



I.3M0). It should 
be noted, however, that its estimated metallicity is 
0.09 (Table ^), being significantly less than as- 



4.3) 



sumed in this section (see 

The ratio C/0 can be > 1 for r ^ 1000 s where the 
abundance of carbon is ~ 0.1 and that of oxygen is 
<^ 0.1. It implies that the carbon-rich ONeMg novae, 
i.e., V1370 Aql and V838 Her, may be explained by 
the models with r ^ 1000 s. Note that carbon tends 
to be ov erproduced in the models with r ^ 1000 s 
(§ 4.2.2 ). This may not, however, change the above 
result with r ^ 1000 s. 

4-2.4- Neon and .sodium 

Neon is the second most abundant metal in the ini- 



tial composition (Table 2^). The abundance of neon 
is not significantly reduced for Tpoak ^ 4 x 10* K, due 
to its rather slow proton capture (Figure |l^). Nev- 
ertheless, the substantial nuclear flow appears in the 
Ne-Na cycle even for Tpoak 2 - 3 x 10* K (Figure |) 
owing to the abundant neon initially present. The 
ratio ^^Ne/'^-'^Ne is clearly correlated with the cooling 
timescale, being small for the shorter r, where the f3~^- 
decay lifetime of ^^Na (~ 32 s) is not negligible. On 
the other hand, the ratio ^'^Ne/^'^Ne is clearly corre- 
lated to the peak temperature, increasing with a rise 



in Tpoak- This is due to the faster proton capture on 
22Ne than on ^ONe. 

The abundance of sodium is ^ 10^"^ for r ^ 100 s, 
due to the steady Ne-Na cycle where ^"Ne is most 



abundant (Figure 11). The isotope ratio is also de- 
termined by their reaction rates as 

22Na/23Na ^ A P^Na(p, a)] /X P^Na(p, 7)] - 10 , 

in the temperature range ~ 2 — 4 x 10* K. On the 
other hand, sodium is abundant 0.01 — 0.1 by 
mass) for r < 100 s, where the /3+-decay lifetimes 
of ^^Mg (~ 6 s) and ^^Mg (~ 16 s) are not neghgible. 
Thus, a part of sodium, which is the decayed product 
of the magnesium isotopes, survives the subsequent 
proton capture. The ratio ^^Na/^'^Na reaches ~ 1, 
owing to the abundant ^^Mg and ^^Mg in the Ne- 
Na region during outbursts. The abundance of ^^Na 
shows a similar trend to that of sodium, clearly cor- 
related to the cooling timescale. This abundance can 
be changed b y the large uncertaint y of the ^^Na(p, 
7)^^Mg rate (Kubono et al. 1994 ; Schmidt et al 



1995| ; |Coc et al. 1995| ; [Kubono et al. 1996|) . However 



it may not be significantly affected for r 100 s, since 
the explosive burning ceases while ^^Mg is abundant. 

The enrichment in neon is characteristic of all the 
observed ONeMg novae. On the other hand, no posi- 
tive detection o f sodium has been reported for recent 
ONeMg novae ( Gehrz et al. 1994 ) , due to lack of use- 
ful lines and, probably, little enrichment in sodium in 
the nova ejecta. An alternative way to check the nu- 
cleosynthesis in the Ne-Na region is to compare with 
the result of the 7-ray line survey of the ^^Na decay 
from a nearby ONeMg nova by CGRO or INTEGRAL 
in the near future. 

4-2.5. Magnesium and aluminum 

Magnesium is one of the abundant elements ini- 
tially present, but rather fragile against proton cap- 
ture. As a result, it is mostly transferred to aluminum 
and silicon via the opened Mg-Al cycle ( Timmcrmann 
et al. 1988 ; Champagne et al. 1988 ). As seen in 
Figure |12|, the abundance of magnesium reaches its 
minimum at r 1000 s, in contrast to carbon (Fig- 
ure |9|). For T ^ 1000 s, it reaches ~ 10~^ due to 
the substantial leakage from the Ne-Na cycle and the 
non-negligible /3+-decay lifetime of ^^Al (~ 10 s). 
Note that the most abundant isotope is always ^^Mg 
due to its slowest proton capture. The isotope ratios 
24]v[g/25j^jg a^j^j 24j^g^26j^jg clearly correlated to 
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the cooling timescale. They are, however, not mono- 
tonic with T but compHcated due to the inflow from 
the Ne-Na cycle and the leakage from the Mg-Al re- 
gion, and the various nuclear paths at high tempera- 
ture (Figure ^ . 

The abundance of aluminum shows a similar trend 
to that of magnesium, correlated to the cooling time- 
scale (Figure 12). The ratio ^^A1/^''A1 is not signifi- 
cantly changed, being close to 

26A1/27A1 = A P^Al(p,7)] /A P^Al(p,7)] ^ 0.1 - 0.5 

in the temperature range ~ 1 — 4 x 10^ K. However, 
the ratio decreases with a reduction in the cooling 
timescale, due to the non-negligible /3+-decay life- 
time of ^^Si (~ 6 s) which is the parent isotope of 
^^Al. Note that, for rather high temperature models 
(Tpeak > 4 X 10^ K), the proton capture on ^^Al is 
faster than its /3+-decay. The subsequent isotope ^^Si 
decays to ^^Mg in ~ 12 s through the isomeric state 
of Al, bypassing its ground state. The double peaks 
in ^^Al (~ 3 X 10~^ by mass) can be seen in Figure [ij. 
The one at lower peak temperatures (~ 1.8 x 10* K) 
is consistent with PSTWS95, STWS98, and JH98, in 
which the abundance of ^^Al decreases with increas- 
ing white dwarf mass. The other peak at higher peak 
temperatures (^ 3 x 10* K) is the consequence of the 
substantial nuclear flux from the Ne-Na region. The 
latter peak, which has not been presented in the pre- 
vious works, is of importance on whether ONeMg no- 
vae can be the significant contributors of the Galactic 
2'^Al. Note that the abundance of ^""Al in the lat- 
ter case does not sub stan tially depend on the initial 
abundance of ^^Mg (§ There are large uncertain- 
ties in the reaction rates of ^^A l(p, 7)^^Si (Wicschcr 
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et 19861 ; |Coc et al. 1995|; [Ihadis et al. 1996]) 



Si(p, 7)^'^P (iHerndl et al. 1995D, ^^ Mgjp, 7)^6 Al 
dCoc et al. 1995t |lhadis et al. 199^, and '^^Aljp 



7)^^Si ( Coc et al. 1995 ; Champagne, Brown, fc Sherr 
19S|3|; Coc et al. 1995). Our trial calculations for a 
few models suggest that these uncertainties change 
the abundance of ^^Al by a factor of ^ 2 — 3. 

The clear dependence of magnesium on the cool- 
ing timescale is useful to constrain (Mwd, -^onv) for 
observed ONeMg novae. The estimated abundance 
of magnesium is ~ 4 x 10~^ — 2 x 10~^ for V693 
CrA, V1370 Aql, and QU Vul (Table |), correspond- 
ing to T < 100 s or r > 10^ s (see Figures ^ and 
p!2[ ). The abundance of aluminum does not signifi- 
cantly vary in the MwD--^^cnv space, being not useful 
to constrain (A/wd, Mcnv)- Nevertheless, the abun- 



dance estimates of aluminum are ~3xl0 '^ — 10 ^for 
V693 CrA, QU Vul, and V351 Pup (Table |), which 
is in good agreement with our results. 

4-2.6. Silicon and phosphorus 

The abundance of silicon reaches ~ 3 x 10~^ for 
Tpcak ^ 2 X 10* K (Figure |l|) via the substantial 
nuclear flux from the Mg-Al region. The abundance 
is only weakly correlated to the cooling timescale. On 
the other hand, the ratios ^SSi/^^Si and ^^Si/^^Si are 
clearly correlated to the cooling timescale, because 
of various competitions between proton capture and 
/J^-decay (Figure ||). 

The abundance of phosphorus ("^^P) reaches ~ 
10-3 - 10-2 for Tpeak ^ 3 X 10* K, due to the faster 
proton capture on ■^'^P than its /J^-decay (Figure p^) . 
Since the Si-P cycle is not closed as seen in Figure 
phosphorus is not significantly destroyed. 

The abundance of silicon in the ejecta of V693 CrA, 
V1370 Aql, and V351 Pup is as smaU as - 2-7 x lO^^, 
corresponding to Tpeak ^ 2 x 10* K. In contrast, that 
in QU Vul (^ 3 — 4 X 10^^) is in agreement with the 
models Tpoak <; 2 x 10* K. The discovery of phospho- 
rus has been reported in the ejected shell of V1974 



Cyg by a near infrared spectroscopy ( [Wagner fc De 



Poy 1996|) . It suggests that V1974 Cyg can be ex- 



plained by the model with a rather high peak temper- 
ature, although an accurate abundance of phosphorus 
is required to constrain (Mwd, Mcnv)- It is also inter- 
esting to note that significantly enhanced phosphorus 
has been detected on the white dwarf in a dwarf nova 
system ( Sion et al. 1997| ) and in the broad line system 
of a QSO (Shields 1996|), which might originate from 



ONeMg novae. 

4.2.7. Sulfur and other heavy elements 

The abundance of sulfur reaches ~ 10-^ for Tpeak ^ 
3 x 10* K, through leakage from the Si-P region (Fig- 
ure HJ). The abundance does not exceed 10-^ in the 
models Mwd ^ 1.15Mq because of the shorter cool- 
ing timescale (Figure This condition is, however, 
highly dependent on the initial composition as will 
be discussed in § For Tpeak ^ 3 x 10* K, the ra- 
'■^^S decrease with a rise in peak 
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tios 32S/33S and ^^S/^ 
temperature, due to the increasing nuclear paths (Fig- 
ure H). For Tpeak ^ 3 X 10* K, these ratios approach 
those determined by their reaction rates. 

At least a half of the observed ONeMg novae, 
V1370 Aql, QU Vul, and V838 Her, are abundant 
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in sulfur in their ejecta (Table H). In addition, the 
sulfur enrichment has been confirmed in the V1974 



wai d et al. 1992 



infrared spectroscopies 


Nood- 


Woodward et al. 1995; 


Wag- 



& DePoy 1996). These novae can be explained 
by the models with such high peak temperatures as 
Tpeak ^ 3 X 10® K. The estimated abundance of sul- 
fur for V1370 Aql is much more abundant than by 
any models in the M^^^-Mcnv space (Figure |l^). It 
should be noted, however, the estimated metallicity 
for V1370 Aql is twice as much as assumed in this 
section (see § 

Heavier elements, from chlorine to calcium, are not 
substantially enhanced for Tp^ak ^ 4 x 10® K (Fig- 
ure . In addition, their enhancement is never seen 
in the models Mwd <^ 1.15Mq due to the shorter 
cooling timescale. Nevertheless, the enrichment in 
chlorine has been reported for the e jecta of V1974 
Cyg by a near infrared spectroscopy (Wagner & De- 
Po f 1996 ). The accurate abundance of chlorine would 
severely constrain (Mwd, -^onv) for V1974 Cyg. 

4.3. Dependence on the initial composition 

So far we have discussed the nucleosynthesis results 
for only one set of the initial composition Xwd = 0.4 
(case B). However, the metallicities of the ejecta for 
V1370 Aql, QU Vul by Saizar et al. (1992), and V838 
Her are significantly deviated from 0.4 (Table In 
addition, the different authors present different metal- 
licity estimates for the same nova events. In partic- 
ular, the discrepancy is serious for QU Vul between 
Saizar et al. (1992) (~ 0.10) and Andrea et al. (1994) 
(~ 0.40). It is, therefore, difficult to judge whether 
the dispersion of the metallicities is real or due to ob- 
servational errors. In the following, we discuss how 
the initial composition influences the nucleosynthesis 
results, comparing the low (X^d =0.1; case A) and 
high {XwD =0.8; case C) metallicity cases. 

As discussed in § |2.l| , the density and temperature 
structures of an envelope are determined uniquely by 
a set of (MwD, -^^cnv) in our model, being indepen- 
dent of its time evolution (but slightly dependent on 
the time variation in mean molecular weight). As a 
result, case C is at most 20 % higher than case A in 
peak temperature for each (AfwDj -^onv) as seen in 
Figure ^ The higher temperature in case C is due 
to the larger mean molecular weight. In contrast, a 
variation in initial composition is crucial for the cool- 
ing timescale (Figure |l|). For Tpoak > 2 x 10® K, 
case C is more than 10 times shorter than case A in 



T. This is a consequence of the higher nuclear energy 
in case C (Figure |l^) due to the abundant nuclear 
fuel. The ejection velocity is also affected by the ini- 
tial composition. As seen in Figure |l^, case C obtains 
significantly higher Vcj than case A in each model. 

A prominent distinction between case A (N0540A) 
and case C (N0540C) can be seen in Figure 17, which 
shows the nuclear flows and the final yields in the 
model (Mwd/M0, M^^v/Mq) = (1.05, 10-^°). In 
N0540A, the nuclear flow extends to sulfur due to the 
longer cooling timescale (~ 23000 s), while that in 
N0540C (r ~ 1800 s) to silicon. The model N0540A 
consumes most of oxygen initially present, in contrast 
to N0540C. 

Figures p^pO| show the abundances of important 
elements and 7-ray emitters in the My^D-Mgnv space 
for case A and case C. These results are explained as 
follows: 

1. The abundance of carbon is still clearly correlated 
to r as in case B (§ 4.2.2|) , reaching its maximum at 
T 1000 s for both cases (Figure |8|). The abundance 
is roughly proportional to Xwd among the models 
with the same cooling timescale. 

2. Magnesium is another element clearly correlated 
to T as in case B (§ 4.2.5 ). In case A, the abundance 
is significantly smaller than in case C, not enhanced 
even for t < 1000 s. This is a consequence of the 
longer r in case A, where the nuclear flow extends to 
heavier elements than magnesium (Figure . 

3. Silicon is also an element showing a correla- 
tion to r in case B, not significantly changed for 
Tpeak ^ 2 X 10® K (§ |T|). This feature holds for 
case C. However, the abundance in case A has a corre- 
lation to Tpoak rather than r, reaching its maximum 
at Tpeak - 2.5 X 10® K (Figure |l|). The depletion 
of silicon in case A for high Tpoak is due to the long 
cooling timescale. 

4. The trend of oxygen abundance significantly dif- 
fers between case A and case C (Figure The 
abundance in case B is correlated to both Tpoak and 
T, being more abundant in the lower Mwd models 
(§ 4.2.3|) . In case C, however, the abundance is not 
significantly changed in the (Mwd, -^^cnv) space, be- 
ing ^ 0.3. On the other hand, that in case A is clearly 
correlated to the peak temperature, significantly de- 
pleted for Tpoak ;> 2.5 X 10® K (Figure ||). 

5. The abundance of sulfur shows a correlation to 
Tpoak in all cases. In case C, however, the abundance 
is < lO^'^ for the models A/wd < 1-15M0 because 
of the shorter r. On the other hand, that in case A 
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reaches 3 x 10~^ at ;> 3 x 10^ K, since r is longer 
and thus the nuclear flow extends to heavier elements. 
6. The radioactive species ^Be, ^^Na, and ^^Al are not 
significantly enhanced in case A because of its longer 
cooling timescale (Figure |2^). On the other hand, 
thebe abundances m case (J show similar trends to 



and V1370 Aql are not discussed in this section, since 
they are not reproduced in our models (§ 4.3). 

shows the models which are in agreement 



Figure 



21 



with the abundance estimates for recent ONeMg no- 



vae, within a factor of three for V693 CrA (Vanland- 



case U m the A^wD^-'Wcnv space (I' lgures 

The estimated metallicity for the ejecta of V1370 
Aql is extremely high (Z ^ 0.85; Table which 
is close to case C. However, the abundance of oxy- 
gen is significantly small (~ 4 — 5 x 10^^), being 
inconsistent with our results 0.1). In addition, 
sulfur in the ejecta is extremely abundant 0.1 by 
mass), which is also in disagreement with our results 
10"^). These features, i.e., the abundances of low 
oxygen and high sulfur, could be explained by lower 
XwD models rather than higher ones (Figure |l^, |lj, 
|8[ and |l^. Thus, this extremely high metallicity 
in this nova ejecta may not be real but due to the 
difficulties in the observational estimations. 



ingham, Starrfield, & Shore 1997; triangles), V351 
Pup (3aizar et al. 1996; asterisks), and V1974 Cyg 
(Austin et al. 1996; stars), and of five for QU Vul 
(Andrea, Drechsel, & Starrfield 1994; circles) and 



^ '"or the QU Vul ejecta. Saizar et al. (1992) gave 



V838 Her ( [Vanlandingham, Starrfield, fc Shore 1997 ; 
squares). The thick symbol for each nova is the best 
model, whose ratio to its observation is shown in Fig- 
ure Interestingly, at least four events (V693 CrA, 
QU Vul, V838 Her, and V1974 Cyg) are well ex- 
plained by the models ~ I.I Mq, which is near the 
lower limit to ONeMg cores (|Nomoto 1984| ). This is 
in contrast to the mass range of 1.25 — 1.35AfQ used 
by PSTWS95 and STWS98, that is near the upper 
bound to ONeMg cores. Table ^ shows the es t imated 
ejecta masses of QU Vul (Taylor et al. 1987 ; [Green- 
house et al. 19881 ^aizar et al. 1992]), V351 Pup 



a r nuch lower metallicity estimate (Z 0.10) corre- 



sp dnding to case A than Andrea et al. (1994). The 



lo w abundance estimates of carbon, oxvgen. and mag- 



nes ium by Saizar et al. (1992) are in good agree- 



Baizar et al. 1996]), V838 Her ( [Woodward et al 
1992[; [Vanlandingham et al. 1996[), and V1974 Cyg 



(Ravelin et al. 1993; Shore et al. 1993; Woodward ei 



ment with our results for Tpeak ^ 2 x 10^ K (Fig- 
ures ^ 1^ and |l^). However, the abundance of sili- 
con (~ 4 x 10~^ by mass) suggests that the nova has 
obtained Tpoak ^ 2 — 3 x 10* K, which is inconsis- 
tent with the above result. Thus, there is no (Mwd, 
Menv) model which explains the abundance estimates 
by Saizar et al. (1992) within reasonable observa- 
tional errors. 

The V838 Her ejecta also shows a rather low metal- 
licity estimate {Z ~ 0.09), which again corresponds to 
case A. The abundance features of the ejected shell, 
i.e., the low oxygen and high sulfur, are well repro- 
duced in our results for Tpoak ~ 2 - 3 x 10* K (Fig- 
ures |l^ and |l^). Hence, the low metallicity for this 
case implies the presence of a real dispersion in metal- 
licity among the observed nova ejecta. 

5. Comparison with observations 



In this section, we discuss which (Af-yvD, M^nv) 
models best match the recent ONeMg nova observa- 
tions from the nucleosynthetic point of view, using the 
resuhs of case B (Xwd = 0.4). For V838 Her, how- 
ever, those of case A (Xwd = 0.1) are used (§ 4.3). 
The abundances for QU Vul by Saizar et al. (1992) 



al. 1997) from observations. These significantly high 
ejecta masses compared with theoretical estimates are 
reasonably explained by our nucleosynthesis results 
if we assume that almost all the envelope is even- 
tually blown off. In addition, for the models with 
Monv ^ 10~'*Mo, the expansion velocities exceed ?;csc 
and obtain v^j ^ 1000 km s^^ (Figures I and |l|), 
which are in good agreement with observations. Note 
that the abundances of carbon and nitrogen by our 
results are also in good agreement with those by ob- 
servations, regardless of their uncertainties (§^). This 
is a consequence that these novae are well explained 
by the models with r <^ 1000 s where the uncertain- 
ties (caused by the depletion of ^"^N) may be small 

(§ H). 



5.1. V693 CrA 

The high oxygen abundance (^ 0.1 — 0.2 by mass) 



in the V693 CrA ejecta ( 


Williams et al. 1985; 


Andrea, 


Drechsel, & Starrfield 1994 




Vanlandingham, Star- 


rfield, & Shore 1997 


) implies 


that it was an event with 



Tpoak < 2 X 10* K or with A/wd 1.15Afo (§ [4.2. 3[ ). 
The low magnesium and high silicon abundances by 
Vanlandingham, Starrfield, & Shore (1997) suggest 
that the cooling timescale was < 1000 s (§ 4.2.5 and 
4.2. 6[). On the other hand, Williams et al. (1985) and 
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Table 5: Ejected Masses of recent ONeMg Novae 





Afej/Af0 


Observations 


QU Vul^ 


8 X 10--* 


Radio emission 


QU Vul2 


> 9 X 10-'' 


Infrared emission 


QU Vul3 


0.2 - 1.5 X 10-'' 


Multiwavelength study 


V351 Pup'i 


1 X 10-^ 


Multiwavelength study 


V838 Her^ 


6.4 - 9 X 10-5 


Infrared emission 


V838 Her^ 


1.8 X lO-'' 


Optical and UV emission 


V1974 Cyg^ 


> 7 X 10-5 


Radio emission 


V1974 Cyg8 


1 - 4 X lO-'' X y-i/2t 


UV emission 


V1974 CygS 


2 - 5 X lO-'' 


Infrared emission 



References: ^Taylor et al. 1987, ^Greenhouse et al. 1988, ^Saizar et al. 1992, "Saizar et al. 1996, ^Woodward et al. 1992, 
''Vanlandingham et al. 1996, '^Pavelin et al. 1993, *Shore et al. 1993, ^Woodward et al. 1997 

ty is the enhancement factor for the helium abundance 



Andrea, Drechsel, & Starrfield (1994) present some- 
what higher magnesium and lower silicon abundances. 
We compare our results with the abundance estimates 
by Vanlandingham, Starrfield, & Shore (1997), since 
others used the overexposed spectrum as pointed out 
by Andrea, Drechsel, & Starrfield (1994). As a re- 
sul t, the model (A/wd/Mq, Menv/Mg) (1.05, lO ''^^ 



10-5-5 _ 10-5), (1.15 _ 1.2, 10"*^ - 10-5-5), and the 
best for (1.25, 10"^) (case B). In such low tempera- 



ture models, magnesium must be abundant (§ 4.2.5), 
though it is not presented in Saizar et al. (1996). 
The above low envelope masses may be due to mass 
accreting at high rate from a giant companion, that is 
also sugges ted by the optical spectral analysis ( ^aizar 



(ca ^e B) is in good agreement with the observation et al. 1996 ). The estimated ejecta mass, 2 x 10 



within a factor of 3 (Figures |21| and^). 

5.2. QU Vul 

The high abundance of sulfur implies that the nova 
obtained such a high temperature as Tpcak > 3 x 10* K 
(Figures || and |l^. Furthermore, the abundance of 
oxygen despite such a high temperature suggests that 
the white dwarf mass was ^ 1.15Mq (§ 1.2.3 ). Our 
results are in agreement with the observational esti- 
mates within a factor of 5 for the models (Mwd /Af© , 
A^env/Afo) = (1.05 - 1.1, 10-3-5 _ 10-3) (case B). 
These high envelope masses are in good agreement 
with the observational estimates of the nova ejecta 
(Table ^ . Note that the both high abundances of oxy- 
gen and sulfur were not explained by previous hydro- 
dynamic studies, with much smaller envelope masses. 

5.3. V351 Pup 

The ejected shell of V351 Pup shows the high oxy- 
gen and low silicon abundances ( Saizar et al. 1996[ ). 
This feature is well explained with the low temper- 
ature models of Tpoak < 2 x 10*^ K (Figures § |l^, 
and p^ . Our results are in good agreement with 
the observational estimates within a factor of 3 for 
the models (Afwo/Afo, M^^^/Mq) = (1.05 - 1.1, 



(Table |5|), implies that this nova occurred in such a 
massive white dwarf as AfwD <; 1.25Af0. 

5.4. V838 Her 

The low oxygen and high sulfur abundances in 
the V838 Her ejecta are the prominent feature in 
the low metallicity models (case A) with Tpeak ~ 
2.5-3x10*^ K (§1^). In addition, the ratios C/N~ 1 
and C/0 ^ 1 suggest that the cooling timescale was 
~ 1000 s (§ [1.2. 2| and |4.2.3D . Thus, the nova may have 
occurred with the low A/wd and the high Afenv (Fig- 
ure |l5| ). Our results are in agreement with the obser- 
vational estimates within a factor of 5 for the model 

(ATwd/Mq, Afenv/Afo) = (1.05, 10-4 - 10-3-5). 

5.5. V1974 Cyg 

Unfortunately, the abundances heavier than neon 
are not presented in Austin et al. (1996), due to 
lack of these lines. The high oxygen abundance sug- 
gests that the peak temperature was ^ 2 x 10* K or 
the white dwarf mass was 1.15A/0 (§ 4.2.3). In 
addition, the ratio C/N ^ 1 implies that the cool- 
ing timescale was ^ 1000 s (§ 1.2.2). Our results 
are in good agreement with the observational esti- 
mates within a factor of 3 for the models (AfwD /Mq , 



14 



Menv/Me) = (1.05, 10-3-5), (1.1, 10-4), and (1.2, 


the solar neighborhood (Starrfield, Truran, & Sparks] 


10- 


5), and the best for (1.1, lO^^.S) (-^^g^, g. no535B, 


197S; 


D'Antona & Matteucci 1991 




Wanajo, Ishi- 


Nli 


)40B, N2050B, and N1045B). They are in reason- 


maru. 


& Kajino 199? ). Recently, some works claim 



able agreement with the estimated mass of the ejecta 



X 10- 



'0 



are also in agreement with the estimates from obser- 
vations ~ 0.75 — l.lAf0 (Paresce et al. 1995; Retter 



that they cannot be the major contributors (Mat- 
(Table [3 ). Their white dwarf masses teucci, D'Antona, & Timmes 1995; JH98) 



Lei Dowitz, 



(Krautter et al 



Ofek 1997D but smaller than 
1996^ 



1.25M0 

The observation shows a 
factor of 2 lower hydrogen abundance than our re- 
sult (Figure |2^). This might be due to the subse- 
quent steady hydrogen burning on the white dwarf as 
pointed out by Krauter et al. (1996). Hayward et al. 
(1996) have derived the neon and magnesium abun- 
dances relative to solar values from a mid-infrared 
observation. If their ratio Ne/Mg ~ 30 is adopted, 
the abundance of magnesium would be 3 x 10-^. 
It favors a relatively high envelope mass model (Fig- 
ure |l^). As a result, N1040B would be the best in this 
case. A recent near infrared measurement has shown 
the presence of the lines of phosphorus and chlorine 
together with sulfur in the V1974 Cyg ejecta (Wagner 
& pePoy 1996| ). This suggests that V1974 Cyg expe- 
rienced Tpoak > 3 X 10^ K. In this case, the higher 
Menv models are also favorable. In addition, the ejec- 
tion velocity in N1040B is ~ 1800 km being good 
agreement with observations (~ 2300 km s-^; Gehrz 
19981), while that in N1045B is 



et bl 
s' 



190 km 



Obviously, further analysis of heavy elements 



is needed to constrain the parameters (A/wd, M^nv) 
for V1974 Cyg. 

6. Production of the radioactive isotopes 

In this section, we discuss the possibilities of de- 
tecting the 7-ray emitters ^Be, ^^Na, and the con- 
tribution to the Galactic ^^Al, based on our nucle- 
osynthesis results in ONeMg novae. Figure |2^ shows 
the total masses of ^Be, ^^Na, and ^^Al produced per 
event for ATwd = 0.4 (case B). As seen in this figure, 
the models A/wd ^ I.IMq with M^nv ^ 10"^^^^ 
which are in good agreement with most of observa- 
tions (§ ^), produce significant amounts of these iso- 
topes. In the rest of this section, all the envelope is 
assumed to be eventually blown off. 

6.1. Gamma-ray emission from '^Be electron 
captures 

Classical novae might be a possible site for the ^Li 
production (the electron-capture products of ^Be) in 



Never- 
theless, the 7-rays (at 478 KeV) from the ''Be elec- 
tron capture would be detectable by CGRO or INTE- 
GRAL ([Hernanz et al. 1996| ; [Jose fc Hernanz 1998 ; 
Gomez- Gomar et al. 1998| ). 



The 7-ray line fiux from the ^Be electron capture 
is estimated as 



F {'Be) 



4 X lO^^cm-^s-^ 

/ x(^Be)Afenv \ 
1^ 5 X j 

g-t/r(^Bo) 



3 kpc 



where d is the distance of the nova system from the 
sun. For the CGRO sensitivity (a few 10-^ cm s-^) 
and a typical distance (d ~ 3 kpc), the mass of ^Be 
~ 5 X lO-^Af0 per event is required to be detected. 
The mass of ''Be per event in our model is over 5-10 



times smaller than required (Figure 23). However, an 
ONeMg nova will be a promising target of ''Be 7-rays 
for INTEGRAL in the near future. Note that CO 
novae may produce about 10 times higher ^Bc than 



ONeM g novae (JII98; Wanajo, Ishimaru, & Kajino 
1998|) . 



6.2. Gamma-ray emission from ^^Na decays 

There has been increasing expectations that an 
ONeMg nova might be the first stellar o bject for 
the detection of the 7-ray emitter 



Truran 1990: Starrfield et al. 1993 



^ ( 


Weiss 


& 


Coc et 


al. 



1995|; PSTWS95; Wanajo 1997a, b; STWS98; JH98; 
Gomcz-Gomar et al. 1998). Nevertheless, no positive 



detection has been reported by COMPTEL on board 
CGRO for the recent ONeMg novae, V351 Pup, V838 



Her, and V1974 Cyg ([yudin et al. 1995) 



As seen in Figure 23, the total mass of ^^Na per 
event (case B) is significantly high in the models 
Menv ^ IO-^Mq. The 7-ray fine flux from the 
decay of ^^Na is estimated as 



F (22Na) 



4 X lO^^cm^^s"^ 

/ x(^^Na)Afenv ^ 
1^ 1 X 10-7A/q j 



3 kpc 
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Table 6: ^^Na Production in ONeMg Novae 





year 


d (kpc) 


XMonv/Mg 


Fo^(cm ''^ s 






V693 CrA 


1981 


11.8 


2.0 X 10"^ 


5.5 X 10- 


4 




QU Vul 


1984 


2.8 


8.9 X 10~^ 


4.4 X 10- 


4 




V351 Pup 


1991 


3.5 


4.0 X 10"" 


1.3 X 10" 


8 


5.5 X 10-5 


V838 Her 


1991 


3.4 


1.9 X 10-8 


6.3 X 10- 


6 


3.3 X 10-5 


V1974 Cyg 


1992 


1.8 


2.5 X 10-8 


3.0 X 10- 


5 


2.3 X 10-5 



"initial flux by this work (best sequence) 
''upper limit by lyudin et al. (1995) 



For the CGRO sensitivity and a typical distance 
(~ 3 kpc), the mass of ^^Na ~ 1O~^M0 per event 
is required. This corresponds to an envelope mass 
of ~ IQ-^M© (Figure ||). Table | shows the dis- 
tances, the masses of ^^Na per event from our results 
(for the best fitted models, see § ||), the expected ini- 
tial 7-ray line fluxes at 1.275 MeV, and the upper 



limits to the COMPTEL observations ( lyudin ct al 
19^ for V693 CrA, QU Vul, V351 Pup, V838 Her, 
and V1974 Cyg. V1370 Aql is omitted here, not ex- 
plained by any (AiwD, Mcnv) models as discussed in 
§ |. V693 CrA and QU Vul may have emitted the 

cm~^s~^. However, 
respectively, 6 x 10~^ 



at present, and will decrease to 
-1 at the launch of INTE- 



7-rays as high as ~ 5 x 10" 
their fluxes have decreased to 
and 1 X 10-5 cm s-^ 
3 X IQ-^ and 5 x 10-^ cm s 
GRAL (~ 2001?). In contrast to the above two no- 
vae, V351 Pup and V838 Her may have yielded much 
lower 7-ray fluxes, which are consistent with the up- 
per limits by COMPTEL. The envelope mass of V351 
Pup might be <, IQ-'^Mq (§ U) so that httle 22Na 
may have been produced. Although V838 Her may 
ha ve o btained a massive envelope such as ~ IQ-'^Mq 
(§ 5.4), the cooling timescale was so long owing to 
the low metallicity (^ 10'* s) that little ^^Na survived. 
The 7-ray flux of the ^^Na decay from V1974 Cyg may 
have been near the sensitivity limit to COM PTEL, 
with the abundance in the best model (§ 5.5) and the 
distance of - 1.8 kpc ( |Chochol et al. 1997| ). Thus, if 
the ejected mass was as high as a few 10"'' M© indeed, 
our model would have produced observable ^^Na (or 
the estimated distance is too short). 

It seems that at least four ONeMg novae (V693 
CrA, QU Vul, V838 Her, and V1974 Cyg) in the past 
twenty years have produced sufficient ^^Na for the 
high sensitivity of INTEGRAL (- 4-5x10-^ cm' s-*). 
The next ONeMg nova in the first decade of the 21st 
century will be a promising candidate for detecting 



the 7-ray emitter ^^Na. 

6.3. Galactic ^^Al production 

Since the discovery by HEAO 3 (Mahoney et al. 
1984), many studies have been carried out to explain 
the presence of ~ 1 — 3Mq of ^^Al in the Galaxy. 
In particular, ONeMg novae have been considered 
to be a promising stellar s it e for the ^^Al produc- 
tion (IWeiss fc Truran 1990|; [Nofar, Shaviv, fc Star- 



rfield 19911 ; ^tarrfield et al. 1993| ; jCoc et al. 1995 
PSTWS95; iKolb & Pohtano 1997; Jose, Hernanz, & 



Coc 1997; STWS98; JH98), as weU as AGB stars 



(Forestini, Paulus, & Arnould 1991), Wolf-Rayet stars 



(Prantzos fc Casse 1986; Meynet et al. 1997|), and 



^ J), 

Type II supernovae ( Walter fc Maeder 1989| ; p^rant- 



zos 1993; Timmes et al. 1995). However, the detailed 
observations of 7- r ay lines at 1.8 Me V by COMPTEL 
( piehl et al. 1994| ; [Diehl et al. 1995| ) have shown that 
the Galactic ^'^Al originates from the youngest stellar 
population associated with the spiral arms and the 



local groups (Prantzos & Diehl 1996; Prantzos 1996) 



This may imply that the major sources of the Galactic 
^^Al are Type II supernovae or Wolf-Rayet stars. 

The mass of ^^Al per event is up to ~ 3 x 1O~'^M0 
in the models M^nv ^ 10""* Af© (Figure |2|). The 
upper limit to the Galactic ^®A1 from ONeMg novae 
is thus estimated as 



3Mo 



\ 40yr 



^)(^ 



/oNcMg \ / X(^^Al)A/oi 



0.25 



3 X 10-'^ 



where i?nova is the nova rate in the Galaxy and 
/oNeMg is the fraction of ONeMg novae. This is in 
good agreement with the estimate from the CGRO 
results. If ONeMg novae are not be the major contrib- 
utors to the Galactic ^^Al, its typical mass per event 
must be somewhat smaller than the above value. 
There are some uncertainties in the Galactic nova 
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rates ([Yungelson, Livio, fc Tutukov 1997 ; Shatter 
19C0; [Hatano et al. 1997| ) and the fraction of O NeMg 
novae dRitter et al. 1991| ; |Livio fc Truran 19941 ). How- 
ever, these uncertainties may be much smaller (a fac- 
tor of ~ 2) than those in the ^^Al yields (about 2 or- 
ders of magnitude as can be seen in Figure |2^) . The 
INTEGRAL survey on the diffuse component of the 
Galactic ^^Al, together with a search of ^^Na from 
an individual ONeMg nova, will impose a severe con- 
straint on the current nova models. 

7. Conclusions 

In this paper we have examined nucleosynthesis 
in ONeMg novae with the wide ranges of the three 
parameters, i.e., the white dwarf mass (Mwd = 
1.05 - 1.35Mq), the envelope mass (Mo„v = 10"'' - 
10~'^Mq), and the initial metallicity (ATwd 0.1 — 

0. 8.. We used a quasi-analytic nova model with a one- 
zone envelope, coupled with an updated nuclear reac- 
tion network code. Our nucleosynthesis results are in 
good agreement with those of previous hydrodynamic 
calculations except for several fragile isotopes. 

We have found that the explosion is more violent 
in a lower Mwd model among those with the same 
peak temperature, due to its smaller gravitational 
potential. There exists the critical cooling timescale 
1000 s), at which the energy generation by the 
/3+-decay of ^''O and ^^O plays a crucial role to the 
envelope expansion. For the models with r < 1000 s, 
the nucleosynthesis results significantly deviate from 
those expected in steady nuclear flows (e.g., the CNO 
and Ne-Na cycle) . These models also obtain high ejec- 
tion velocities (> 1000 km s""'^), which are in good 
agreement with recent observations. 

There are a couple of characteristic trends for 
the abundances in the Mwd-Mohv space as follows 
(case B): 

1. The abundances of oxygen, neon, phosphorus, and 
sulfur (and ''Be) are clearly correlated to the peak 
temperatures, although those of oxygen and sulfur are 
also dependent of the cooling timescales. The abun- 
dance of oxygen is always abundant in the models 
■Mwd <^ 1.15M0. The heavier elements than sulfur 
show no significant enrichment in the models with 

Tpeak ^ 4 X 108 K. 

2. The abundances of carbon, fluorine, sodium, and 
magnesium (and ^^Na, ^^Al) are clearly correlated to 
the cooling timescales. The abundance of ^^Na is sig- 
nificantly high in the models with r < 100 s. On the 



other hand, that of ^^Al shows double peaks in the 
MwD-Mcnv space. 

3. The abundances of nitrogen, aluminum, and sili- 
con are not significantly changed in the MwD-Afenv 
space, although those are weakly dependent of the 
cooling timescales. 

The initial metallicity A'wd, as well as Mwd and 
Mcnv, is a crucial parameter to the nucleosynthesis 
results. For smaller ATwd, the explosion is less vio- 
lent and thus the cooling timescale is longer, because 
of the smaller nuclear fuel. As a result, the models 
with low ATwd (case A) produce more sulfur but less 
oxygen than those with high ATwd (cases A and B). 
The former case is unfavorable for the production of 
^Be, 22Na, and ^SAl. 

Comparison of our nucleosynthesis results with ob- 
servational abundance estimates enables us to con- 
strain the model parameters (Mwd, Mcnv) for the 
observed ONeMg novae. We have found that the 
white dwarf masses of at least four of the observed 
six ONeMg novae are as low as ~ I.IMq. This is 
significantly smaller than the prediction of Mwd ~ 
1.25 — 1.35Mq obtained by previous hydrodynamic 
studies. On the other hand, our results suggest that 
their envelope masses were ^ 1O"^M0 which are 
consistent with the observational estimates of their 
ejected masses. In addition, the observed fast ejec- 
tion velocities for these novae (^ 1000 km s"^) are 
also obtained in those models. There remains a dis- 
crepancy between these high ejected masses and those 
estimated by previous hydrodynamic studies. How- 
ever, a low mass white dwarf (Mwd — I.IMq) may 
be able to accumulate such a massive envelope with a 
small mass accretion rate and a low surface tempera- 
ture (Starrfield et al. 1998). 

Our results also show that the models Mwd — 
I.IMq with Menv ^ 10~^Mq, which are the possible 
explanations for most of the observed ONeMg novae, 
produce significant amounts of ^Be, ^^Na, and ^^Al: 

1. The 7-ray line flux from the ''Be electron capture 
is too weak to be detected with the CGRO sensitiv- 
ity for the typical distance from the sun. However, a 
nearby ONeMg nova could emit the 7-rays detectable 
by INTEGRAL in the near future. 

2. The mass of ^^Na per event is significantly high in 
the models Mo„v ^ 10~'^ Mq. V1974 Cyg may have 
produced an interesting amount of ^^Na which is near 
the upper limit to the COMPTEL sensitivity. Fur- 
thermore, we suggest that at least four ONeMg novae 
in the past twenty years have produced enough ^^Na 
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to the INTEGRAL sensitivity. The next ONeMg nova 
will be a promising target for the detection of the 7- 
ray emitter ^^Na. 

3. The mass of ^*'A1 per event is also significantly 
high in the models M^nv ^ 10"'' M©. The mass of 
the Galactic ^^Al which originates from ONeMg no- 
vae is estimated to be <; 3Mq, being consistent to 
the COMPTEL result. They may not be, however, 
major contributors according to the 7-ray survey at 
1.8 MeV by COMPTEL. The 7-ray line survey by 
INTEGRAL will significantly constrain the ranges of 
(MwD, Menv) for ONeMg novae. 

We should emphasize that hydrodynamic investi- 
gations including multi-dimensional calculations, es- 
pecially with a massive envelope, are necessary to 
prove our conclusions with the one-zone approxima- 
tion. There are also other observables besides the 
abundances which cannot be dealt with in this study 
(e.g., the surface luminosity and the ejecta masses). 
Nevertheless, our results afford some new perspectives 
on the future nova modelings. The future INTE- 
GRAL survey for the 7-ray emitters, together with 
abundance analyses by ultraviolet, optical, and near 
infrared spectroscopies, will also impose a severe con- 
straint on the current nova models. 

We would like to acknowledge useful discussions 
with T. Kajino, S. Kubono, I. Hachisu, and J. W. 
Truran. We would like to express sincere appreciation 
to F. -K. Thielemann for providing the data of nuclear 
reaction rates. This work has been supported in part 
by the grant-in- Aid for Scientific Research (05242102, 
06233101) and COE research (07CE2002) of the Min- 
istry of Education, Science, and Culture in Japan, 
and from Japan Society for Promotion of Science. 



REFERENCES 

Anders, E., & Grevesse, N. 1989, Geochim. Cos- 
mochim. Acta, 53, 197 

Andrea, J., Drechsel, H., & Starrfield, S. 1994, A&A, 291, 
869 



Austin, S. J., Wagner, R. M., Starrfield, S., Shore, S. N., 
Sonneborn, G., & Bertram, R. 1996, AJ, 111, 869 

Boffin, H. M. J., Paulus, G., Arnould, M., & Mowlavi, N. 

1993, A&A, 279, 173 

Caughlaln, G. R., & Fowler, W. A. 1988, Atomic Data 
Nucl. Data Tables, 40 

Gehrz, R. G., Truran, J. W., Williams, R. E., & Starrfield, 

S. 1998, PASP, 110, 3 

Champagne, A. E., Cclla, C. H., Kouzcs, R. T., Lowry, M. 
M., Magnus, P.V., Smith, M. S., & Mao, Z. Q. 1988, 
Nucl. Phys. A, 487, 433 

Champagne, A. E., Brown, B. A., & Sherr, R. 1993, 
Nucl. Phys. A, 556, 123 

Chochol, D., Grygar, J., Pribulla, T., Komzik, R., Hric, 

L., & Elkin, V. 1997, A&A, 318, 908 

Coc, A., Mochkovitch, R., Oberto, Y., Thibaud, J. P., & 
Vangioni-Flam, E. 1995, A&A, 299, 479 

D'Antona, F. & Matteucci, F. 1991, A&A, 248, 62 

Dichl, R,., Dnpraz, C, Bennet, K., Bloomcn, H., dcDoer, 
H., Hcrmson, W., Lichti, G. G., McConnell, M., Morris, 
D., Ryan, J., Schonfeldcr, V., Stcinle, H., Strong, A. 
W., Swanenburg, B. N., Varendorff, M., & Winkler, C. 

1994, ApJS, 92, 429 

Diehl, R., Dupraz, C., Bennet, K., Bloemen, H., Hermsen, 
W., Knodlseder, J., Lichti, G., Morris, D., Ryan, J., 
Schonfelder, V., Steinle, H., Strong, A. W., Swanen- 
burg, B. N., Varendorff, M., & Winkler, C. 1995, A&A, 

298, 445 

Forestini, M., Paulus, G., & Arnould, M. 1991, A&A, 252, 
597 

Fujimoto, M. Y. 1982a, ApJ, 257, 752 
. 1982b, ApJ, 257, 767 

Gehrz, R. D., Woodward, C. E., Greenhouse, M. A., Star- 
rfield, S., Wooden, D. H., Witteborn, F. C, Sandford, 
S. A., AUamandola, L. J., & Bregman, J. D. 1994, ApJ, 
421, 762 

Glasner, S. A., Livne, E., & Truran, J. W. 1997, ApJ, 475, 

754 

Gomez-Gomar, J., Hernanz, M., Jose, J., & Isern, J 1998, 
MNRAS, 296, 913 

Greenhouse, M. A., Grasdalen, G. L., Hayward, T. L., 
Gehrz, R., D., & Jones, T. J. 1988, AJ, 95, 172 

Hamada, T. & Salpctcr, E. E. 1961, ApJ, 134, 683 

Hashimoto, M., Iwamoto, K., & Nomoto, K. 1993, ApJ, 
414, 105 

Hatano, K., Branch, D., Fisher, A. & Starrfield, S. 1997, 
MNRAS, 290, 113 



18 



Hayward, T. L., Saizar, P., Gehrz, R. D., Benjamin, R. 
A, Mason, C. G, Houck, J. R., Miles, J. W., Gull, G. 
E., & Schoenwald, J. 1996, ApJ, 469, 854 

Hernanz, M., Jose, J., Coc, A., & Isern, J. 1996, ApJ, 465, 
27 

Herndl, H., Gorres, J., Wiescher, M., Brown, B. A., & 
Van Wormer 1995, Phys. Rev. C, 52, 1078 

Hillebrandt, W. & Thielemann, F. -K. 1982, ApJ, 255, 
617 

Iben, I. Jr., & Tutukov, A. V. 1985, ApJS, 58, 661 

Iben, I. Jr., Pujimoto, M. Y., & MacDonald, J. 1991, ApJ, 

375, 27 

Ibcu. I. Jr., & Tutukov, A. V. 1996, ApJS, 105, 145 

Icliiruaru, S. & Kitamura, H. 1994, in Proc. Int. Semi, on 
Elementary Processes in Dense Plasmas, ed. S. Ichi- 
maru & S. Ogata (Addison- Wesley) , 113 

Iglesias, C. A., & Rogers, F. J. 1993, ApJ, 412, 752 

Iliadis, C., Buchmann, L., Endt, P. M., Herndl, H., & 

Wiescher, M. 1996, Phys. Rev. C, 53, 475 

lyudin, A. F., Bennett, K., Blocmcn, H., Diehl, R., 
Hermsen, W., Licliti, G., Morris, D, Ryan, J., 
Schondelder, V., Steinlo, H., Strong, A., Varendorff, 
M., & Winkler, C. 1995, A&A, 300, 422 

Jose, J., Hernanz, M., & Coc, A. 1997, ApJ, 479, 55 

Jose, J., & Hernanz, M. 1998, ApJ, 494, 680 

Kercek, A., Hillebrandt, W., & Truran, J. W. 1998a, 

A&A, 337, 379 

-. 1998b, A&A, submitted 

Kolb, U. & Politano, M. 1997, A&A, 319, 909 

Kovetz, A. & Prialnik, D. 1997, ApJ, 477, 356 

Krauttcr, J., Ogclman, H., Starrficld, S., Wichmann, R., 
& Pfcffcrmami, E. 1996, ApJ, 456, 788 

Kubono, S., Yun, C. C., Boyd, R. N., Buchmann, L. 
R., Fuchi,Y., Hosaka, M., Ikeda, N., Jiang, C. L., 
Katayama, I., Kawashima, H., Miyatake, H., Niizeki, 
T., Nomura, T., Odahara, A., Ohura, M., H. Orihara, 
H., Rolfs, C., Shimoda, T., Tajima, Y., Tanaka, M. H., 
& Toyokawa, H. 1994, Z. Phys. A, 348, 59 

Kubono, S., Hosaka, M., Strasser, P., Guimaraes, V., 

Jeong, S. C., Katayama, I., Miyachi, T., Nomura, 
T., Tanaka, M. H., Fuchi, Y., Kawashima, H., Kato, 
S., Yun, C. C., Orihara, H., Niizeki, T., Miyatake, 
H., Shimoda, T., Kajino, T., and Wanajo, S. 1997, 
Nucl. Phys. A, 621, 195 

Kutter, G. S., & Sparks, W. M. 1987, ApJ, 321, 386 

Livio, M., & Truran, J. W. 1994, ApJ, 425, 797 

MacDonald, J. 1983, ApJ, 267, 732 



Mahoncy, W. A., Ling, J. C., Wheaton, W. A., & Jacob- 
son, A. S. 1984, ApJ, 286, 578 

Matteucci, F., D'Antona, F., & Timmes, F. X. 1995, 
A&A, 303, 460 

Meynet, G., Arnould, M., Prantzos, N., & Paulus, G. 
1997, A&A, 320, 460 

Nofar, I., Shaviv, G., & Starrfield, S. 1991, ApJ, 369, 440 

Nomoto, K. 1984, ApJ, 277, 291 

-. 1987, ApJ, 322, 206 

Nomoto, K., & Hashimoto, M. 1988, Phys. Rep., 163, 13 

Politano, M., Starrfield, S., Truran, J. W., Weiss, A., & 
Sparks, W. M. 1995, ApJ, 448, 807 

Paresce, F., Livio, M., Hack, W., & Korista, K. 1995, 

A&A, 299, 823 

Ravelin, P. E., Davis, R. J., Morrison, L. V., Bode, M. F., 
& Ivison, R. J. 1993, Nature, 363, 424 

Prantzos, N., & Cassc, M. 1986, ApJ, 307, 324 

Prantzos, N. 1993, ApJ, 405, L55 

Prantzos, N., & Diehl, R. 1996, Phys. Rep., 267, 1 

Prantzos, N. 1996, A&AS, 120, 303 

Prialnik, D., & Kovetz, A. 1984, ApJ, 281, 367 

Prialnik, D., & Kovetz, A. 1995, ApJ, 445, 789 

Retter, A., Lcibowitz, E. M., & Ofek, E. O. 1997, MN- 
RAS, 286, 745 

Ritossa, C., Garcia, E., -B., & Iben, I. 1996, ApJ, 460, 
489 

Ritter, H., Politano, M., Livio, M., & Webbink, R. F. 
1991, ApJ, 376, 177 

Saizar, P., Starrfield, S., Ferland, G. J., Wagner, R. M., 
Truran, J. W., Kenyon, S. J., Sparks, W. M., Williams, 
R. E., & Stryker, L. L. 1992, ApJ, 398, 651 

Saizar, P., Pachoulakis, I., Shore, S. N., Starrfield, S., 
Williams, R. E., Rothschild, E., & Sonneborn, G. 1996, 
MNRAS, 279, 280 

Schmidt, S., Rolfs, C., Schulte, W. H., Trautvetter, H. P., 

Kavanagli, R. W., Hategan, C., Faber, S., Valnion, B. 
D., & Graw, G. 1995, Nucl. Phys. A, 591, 227 

Shafter, A. W. 1997, ApJ, 487, 226 

Shields, G. A. 1996, ApJ, 461, 9 

Shore, S. N., Sonneborn, G., Starrfield, S., Gonzalez- 
Riestra, R., & Akc, T. B. 1993, AJ, 106, 2408 

Sion, E. M., Cheng, F. H., Sparks, W. M., Szkody, P., 
Huang, M., & Hubeny, I. 1997, ApJ, 480, 17 

Snijders, M. A. J., Batt, T. J., Roche, P. F., Seaton, M. 
J., Morton, D. C, Spoelstra, T. A. T., & Blades, J. C. 
1987, MNRAS, 228, 329 



19 



Starrficld, S., Truran, J. W., & Sparks, W. M. 1978, ApJ, 

226, 186 

Starrfield, S., Sparks, W. M., & Shaviv, G. 1988, ApJ, 
325, 35 

Starrficld, S., Truran, J. W., Politano, M., Sparks, W. M., 
Nofar, I., & Shaviv, G. 1993, Pliys. Rep., 227, 223 

Starrfield, S., Truran, J. W., Wiescher, M. C., Sparks, W. 

M. 1998, MNRAS, 296, 502 

Sugirnoto, D., & Fujimoto, M. Y. 1978, PASJ, 30, 467 

Taylor, A. R., Seaquist, E. R., Hollis, J. M., & Pottasch, 
S. R. 1987, A&A, 183, 38 

Timmcs, F. X., Woosley, S. E., Hartmann, D. H., & Hoff- 
man, R. D. 1995, ApJ, 449, 204 

Timmernamm, R., Becker, H. W., Rolfs, C., Schroder, U., 
& Trautvetter, H. P. 1988, Nucl. Phys. A, 477, 105 

Thielemann, F.-K., Kratz, K.-L., Pfeiffer, B., Rauschcr, 
T., van Wormer, L., & Wiescher, M. C. 1994, 
Nucl. Phys. A, 570, 329 

Thielemann, F.-K. 1995, private communication 

Truran, J. W., Starrfield, S., Strittmatter, P. A., Wyatt, 
S. P., & Sparks, W. M. 1977, ApJ, 211, 539 

Truran, J. W. 1982, Essays in Nuclear Astrophysics, ed. 
C. A. Barnes, D. D. Clayton, & D. N. Schramm (Cam- 
bridge: Cambridge University Press), p. 467 

Truran, J. W., Thielemann, F.-K., & Arnould, M. 1987, 

Nuclear Astrophysics, od. Hillobrandt (Springer), p. 91 

Vanlandingham, K. M., Starrfield, S., Wagner, R. M., 
Shore, S. N., & Sonneborn, G. 1996, MNRAS, 282, 563 

Vanlandingham, K. M., Starrfield, S., & Shore, S. N. 1997, 

MNRAS, 290, 87 

Van Wormer, L., Gorres, J., Iliadis, C, Wiescher, M., & 
Thielemann, F.-K. 1994, ApJ, 432, 326 

Wagner, R. M., & DePoy, D. L. 1996, A&A, 218, 123 

Walter, R., & Maeder, A. 1989, A&A, 218, 123 

Wanajo, S., Nomoto, K., Hashimoto, M., Kajino, T., & 
Kubono, S. 1997, Nucl. Phys. A, 616, 91 

Wanajo, S., Nomoto, K., Truran, J. W., & Hashimoto, M. 
1997, Nucl. Phys. A, 621, 499 

Wanajo, S., Ishimaru, Y., & Kajino, T. 1998, in Origin of 
Matter And Evolution of Galaxies 97, ed. S. Kubono, 
T. Kajino, K. Nomoto, & I. Tanihata (World Scien- 
tific), 262 

Ward, R. A., & Fowler, W. A. 1980, ApJ, 238, 266 

Weiss, A., & Truran, J. W. 1990, A&A, 238, 178 

Wiescher, M., Gorres, J., Thielemann, F.-K., & Ritter, H. 
1986, A&A, 160, 56 



WiUiams, R. E., Noy, E. P., Sparks, W. M., Starrfield, S. 
G., Wyckoff, S., & Truran, J. W. 1985, MNRAS, 212, 
753 

Woodward, C. E., Gchrz, R. D., Jones, T. J., & Lawrence, 

G. F. 1992, ApJ, 384, 41 

Woodward, C. E., Greenhouse, M. A., Gehrz, R. D., 
Pendleton, Y. J., Joyce, R. R., Van Buren, D., Fis- 
cher, J., Jennerjohn, N. J., & Kaminski, C. D. 1995, 
ApJ, 438, 921 

Woodward, C. E., Gehrz, R. D., Jones, T. J., Lawrence, 
G. F., & Skrutskie, M. F. 1997, ApJ, 477, 817 

Yungelson, L., Livio, M., & Tutukov, A. 1997, ApJ, 481, 
127 



This 2-column preprint was prepared with the AAS lAT^^X 
macros v4.0. 



20 



Fig. 1. The (AfwD, ^'^cnv) scquciiiccs at which our 
numerical calculations have been carried out. The 
dots denote this work, and squares, triangles, and 
stars are, respectively, from Politano et al. (1995), 
Starrfield et al. (1998), and Jose and Hernanz (1998). 

Fig. 2. — Contours of the proper values for the pres- 
sure (Ph/fh) and the density (pb/^/b) in the loga- 
rithmic scale in the MwD-Afenv space. 

Fig. 3. The Af i? relations for various white dwarfs. 
The solid line is for O : Ne = 5 : 3 (this work), 
dotted for the completely degenerate electron gas by 
Chandrasekhar's method (Yc = 0.5), broken and dot- 
dash for carbon and magnesium (Hamada & Salpeter 
1961). The dots on the lines for O -|- Ne, carbon, 
and magnesium denote at which neutronization oc- 
curs. The triangles are taken from PSTWS95 and 
STWS98. 

Fig. 4. — The ratios of our nucleosynthesis results to 
STWS98 (sequence 6). The dots and triangles denote 
isotopes and elements, respectively. 

Fig. 5. — The nucleosynthesis results for several 
{M^^b/Mq, Menv/-^©) sequences in the N-Z plane. 
The size of a circle indicates the yield at the final 
stage, and the length of an arrow the net nuclear flow 
in the logarithmic scale. The initial compositions are 
shown by dotted circles. 

Fig. 6. — Contours of the peak temperatures at the 

base, the cooling timescales, the energy generation 
rates per unit mass, and the ejection velocities in the 
MwD-Afenv space (case B). 

Fig. 7. — The time variations of the temperature at 
the base and the nuclear energy generation rate per 
unit mass for (Mwd/Mq, M^^v/Mq) = (1.15, IQ-^-O) 
and (1.35, IQ-^ '^). 

Fig. 8. — Contours of the abundances of ^Be and ^^B 
in the logarithmic scale (case B). 

Fig. 9. — Same as Figure 8, but for carbon, nitrogen, 
and their isotope ratios. 

Fig. 10. — Same as Figure 8, but for oxygen, its iso- 
tope ratios, and fluorine. 

Fig. 11. — Same as Figure 8, but for neon, sodium, 
their isotope ratios, and ^^Na. 



Fig. 12. Same as Figure 8, but for magnesium, 
aluminmn, their isotope ratios, and ^^Al. 

Fig. 13. Same as Figure 8, but for silicon, its isotope 
ratios, and phosphorus. 

Fig. 14. — Same as Figure 8, but for sulfur and its 
isotope ratios, and the sum of chlorine, argon, potas- 
sium, and calcium. 

Fig. 15. — Contours of the peak temperatures at the 
base and the cooling timescales in the MwD-Afenv 
space for case A and C. 

Fig. 16. — Same as Figure 15, but for the energy gen- 
eration rates per unit mass and the ejection velocities. 

Fig. 17. — The nucleosynthesis results for (AfwD /Mq , 
Afenv/Me) = (1.05, 10-''-°) in the N-Z plane for 
case A and C. 

Fig. 18. — Contours of the abundances of carbon and 
oxygen in the logarithmic scale for case A and C. 

Fig. 19. — Same as Figure 18, but for magnesium, 
silicon, and sulfur. 

Fig. 20. — Same as Figure 18, but for '^Be, ^^Na, and 

26 Al. 

Fig. 21. — The (AfwD, Afonv) sequences which are 
in agreement with recent ONeMg novae, within the 
factor of three for V693 CrA (triangles), V351 Pup 
(asterisks), and V1974 Cyg (stars), and of five for 
QU Vul (circles) and V838 Her (squares). The thick 
signs are the best sequences in our results. 

Fig. 22. — The ratios of our results to observational 
abundance estimates. The symbols are the same as 
Figure 20. 

Fig. 23. — Contours of the masses of ^Be, ^^Na, and 
26 Al per event in the logarithmic scale (case B). 
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